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Preface 



Aeronomy , a part of atmospheric science, is defined as the study of the physics 
and chemistry of the upper atmosphere of a planet. The name was coined 
about 50 years ago by a pioneer of this discipline, the British geophysicist 
Sydney Chapman, to differentiate it from meteorology , which is concerned 
with the lower atmosphere where the prevailing weather phenomena can be 
described by bulk properties of a gas. In contrast to this, aeronomy, dealing 
with the tenuous upper part of the atmosphere, must apply a gas kinetic 
approach, describing processes such as molecular diffusion and heat conduc- 
tion as well as photochemical processes such as dissociation and ionization of 
gases. 

Sydney Chapman himself is credited with such fundamental contributions 
to aeronomy as the formation of stratospheric ozone and the theory of iono- 
spheric layer formation, which carries his name, as well as with the highly 
useful concept of atmospheric scale height. Beyond photochemistry, aeron- 
omy is also concerned with atmospheric escape by thermal and non-thermal 
processes that are playing an important role in the evolution of planetary 
atmospheres. 

The present book is not intended as a traditional textbook but rather 
as a compendium of the most important physical and chemical processes 
occurring in planetary atmospheres. Most of the material is taken in slightly 
modified or corrected form from the monograph by one of us (SJB), Physics of 
Planetary Ionospheres , published by Springer Verlag thirty years ago (1973) 
as Vol. 6 of the series “Physics and Chemistry of Space” , which had a similar 
purpose: to address the physical and chemical processes that lead to the 
formation of ionized layers, as well as to the heating of the thermosphere as 
the result of the interaction of ionizing radiation with a neutral atmosphere. 
The references cited in the text therefore are a selective mix of “old” and new 
contributions with a separate bibliography containing comprehensive reviews 
of topics covered by this monograph. 

Just as the “old” book was intended to concentrate on the underlying 
physical processes for explaining observations then expected as the result of 
the early exploration by space probes of the planets and moons in our Solar 
System, the new version, called Planetary Aeronomy , aims at providing in 
succinct form also the physical background to be applied to information about 
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the atmospheric environments and evolution of planets outside our Solar 
System, the so-called exoplanets, whose atmospheres will become observable 
at an increasing rate during the next few decades by space missions like ESA’s 
Darwin and NASA’s Terrestrial Planet Finder (TPF) (launch foreseen in the 
second decade of this century) and their precursors like CoRoT (Convection 
Rotation and planetary Transits; CNES, launch 2006), Kepler (NASA, launch 
2007), Eddington (NASA), etc. that are currently under development. 

Our book does not contain sophisticated models , some already in the liter- 
ature, but rather limits itself to their physical background. It is in consonance 
with an adage attributed to Aristotle: It is the mark of an instructed mind to 
rest satisfied with the degree of precision which the nature of the subject per- 
mits and not to seek an exactness where only an approximation of the truth 
is possible! 

Graz, February 2004 S.J. Bauer and H. Lammer 
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Introduction 



For almost 50 years investigations of the Earth’s upper atmosphere have 
provided comprehensive knowledge about the subdiscipline of atmospheric 
science, now generally known as aeronomy. Since the advent of the space age 
observations made by spacecraft have made some of the traditional sciences 
into a “space science”. This is particularly true of the possibilities created 
for astronomical observations outside our Earth’s atmosphere with complex 
astronomical observatories, which have led to an unprecedented increase of 
information about stellar radiation sources beyond our own Sun. 

While “visible” radiation from a stellar source usually penetrates the at- 
mosphere of a planet down to its surface (unless blocked by clouds and dust 
particles) , the “short wave” radiation of high photon energy in the extreme 
ultraviolet and X-ray (XUV) part of the spectrum is generally absorbed in 
the more tenuous part of the upper atmosphere by its constituents, leading to 
dissociation of molecular species, ionization and heating. Thus, properties of 
this radiation (including its long term evolution) and the neutral atmosphere 
with which it interacts, as well as the results of this interaction, i.e., heating 
leading to its characteristic temperature, the formation of ionized layers and 
the attendant photochemistry as well as transport processes of the neutral 
and ionized species, have to be considered in the treatment of aeronomical 
processes. 

As the result of the low particle concentrations reached at the high alti- 
tudes on account of the exponential decay of pressure and density, collisions 
among gas particles become increasingly less important, so that a gas kinetic 
theory approach has to be used. In the limit the upper part of a gravitation- 
ally bound atmosphere, the “exosphere” becomes essentially collisionless, so 
that particles arising from below with sufficient velocities can even escape and 
this leads to a loss of atmospheric constituents to space due to their thermal 
motion. Beyond this, escape can also occur via photochemical processes and 
plasma processes (i.e., interaction of stellar winds with ionospheric plasma, 
when there is no protective planetary magnetic field). 

In the following chapters we discuss in a succinct fashion the relevant 
properties and processes, keeping in mind the more general case of planets 
(and their satellites) even beyond our own Solar System, the so-called ex- 
oplanets, and their stellar energy sources. Although most of our knowledge 
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about aeronomical processes has been gained from the study of our own Solar 
System, they obviously can be applied to extrasolar planetary systems. 

In spite of the broad general adherence to the “International System” of 
units (SI), some scientific disciplines (including our own) are still using units 
based on the old cgs system (particularly the length scale centimeters instead 
of meters) even in the recent published international literature. In this book 
we therefore adhere to the same custom, although providing the appropriate 
conversions in the text and the appendix. In addition commonly used units 
in astronomy and aeronomy as well as numerical values for planetary bodies 
are also given in the appendix. 




I Radiation and Particle Environments 
from Stellar and Cosmic Sources 



The principal ionizing radiation responsible for the heating of the upper plan- 
etary atmospheres and the formation of planetary ionospheres from stellar 
or solar sources are extreme ultraviolet and X-rays (XUV) and corpuscular 
radiation, i.e., galactic cosmic rays and energetic particles having their origin 
in stellar systems (cosmic ray protons, auroral particles, stellar wind). 

The atmosphere of a planet can be divided into a number of regions. It 
is quite natural that the terminology which has evolved for the terrestrial 
atmosphere should also be applied to other planetary atmospheres. The ter- 
minology which is generally accepted is based on the vertical distribution of 
temperature and composition, i.e., the parameters necessary for describing 
an atmosphere for which the hydrostatic condition can be used. 

In terms of the temperature variation, a planetary atmosphere can be 
divided as follows (see Fig. 1.1). The lowermost part of the atmosphere is 
the troposphere , where the primary heat source is the planetary surface and 
heat is convected by turbulent motion. This leads to a convective or adiabatic 
temperature distribution. The vertical temperature gradient dT/dz , or adi- 
abatic lapse rate is given by dT/dz — — g/c p , where g is the acceleration of 
gravity and c p is the specific heat at constant pressure. It is therefore depen- 
dent on the planet’s acceleration of gravity and the atmospheric composition. 
For Earth, the observed adiabatic lapse rate is about —6.5 K/km due to the 
presence of water vapor and large-scale planetary motions. 

The troposphere terminates at the tropopause , the level at which the tem- 
perature decrease with an adiabatic lapse rate ceases. Above the tropopause 
the temperature distribution is governed by radiative rather than convective 
processes and the temperature decreases much more slowly (\dT/dz\ < g/c p ) 
or becomes essentially constant ( dT/dz ~ 0). For Earth the tropopause oc- 
curs at an altitude of about 13 km. Above the actual tropopause is the strato- 
sphere. This region was originally thought to be isothermal throughout. In 
the terrestrial stratosphere the temperature, after being initially constant, 
increases with altitude due to the ultraviolet (UV) absorption by ozone (O 3 ), 
reaching a temperature maximum at the stratopause (about 50 km). Above 
this level the mesosphere begins, where dT/dz < 0, reaching a temperature 
minimum at the mesopause (Earth: about 85 km), due to the presence of CO 2 
and H 2 0 which provide a heat sink by radiating in the infrared. Above the 
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Fig. 1.1. Nomenclature of atmospheric regions based on terrestrial analogies. 



mesopause, XUV radiation is absorbed and in part is used for heating, lead- 
ing to a positive temperature gradient, dT/dz > 0. This region is called the 
thermosphere . In the lower thermosphere convection is the principal process 
of heat transport, while in the upper thermosphere heat is transported by 
conduction, leading to an isothermal region (T = const.). The start of the 
isothermal region is termed the thermopause. In this region, the mean free 
path of the gas becomes large and collisions become negligible, so that light 
atmospheric constituents, whose velocity exceeds the gravitational escape ve- 
locity, Uoo, can thermally escape from the planet. 

This region is called the exosphere ; the isothermal temperature above the 
thermopause is therefore also referred to as the exospheric temperature. The 
exosphere begins at an altitude which was originally called the critical level , 
but which is now generally referred to as the exobase. This level is defined 
by the condition that the mean free path is equal to the local scale height, 
i.e., the logarithmic decrement of pressure with altitude. The exobase may 
also be called the baropause , since the entire atmosphere below that level is 
also referred to as the barosphere , i.e., the region where the barometric law 
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holds perfectly. In the exosphere the velocity distribution is non-Maxwellian 
due to the escape of the high velocity particles and the density does not 
strictly follow the barometric formula, but has to be derived by considering 
the individual ballistic components of the atmospheric gas. 

The troposphere and stratosphere together are also referred to in a less 
precise way as the lower atmosphere , while the regions above the stratosphere 
are called the upper atmosphere. It should be noted that without a strato- 
spheric heat source, such as the absorbtion of O 3 in the case of Earth, a 
planetary atmosphere may not possess a stratopause and the region above 
the tropopause may be called either the stratosphere or the mesosphere. 
The lower atmosphere is the principal domain of meteorology whereas the 
upper atmosphere is that of aeronomy. However, it would be a mistake to 
rigidly compartmentalize the atmosphere, since interactions between the var- 
ious regions do occur. For the discussion of planetary aeronomic processes, 
the emphasis must be, however, on the upper atmosphere. 

The atmosphere can also be divided in terms of its composition into the 
homosphere and the heterosphere. These terms are somewhat less frequently 
used than the nomenclature based on the temperature distribution. In the 
homosphere the composition is uniform (Earth: 78% N 2 , 21% O 2 , 1% Ar 
plus minor constituents). This region is characterized by turbulent mixing. 
In the heterosphere, the composition varies due to dissociation of molecu- 
lar constituents and as the result of diffusive separation. The level where 
diffusion rather than turbulent mixing becomes the controlling process is 
generally called the homopause (formerly also called the turbopause). More 
precisely, the homopause can be defined as the level where the eddy diffu- 
sion (mixing) coefficient is equal to the molecular diffusion coefficient. Since 
these coefficients are different for different constituents, they will have differ- 
ent homopause levels. In the terrestrial atmosphere, homopause levels occur 
at altitudes of about 100 km. The homopause concentrations of atmospheric 
constituents represents important boundary conditions for their distribution 
in the upper atmosphere and the ionosphere. 

The ionosphere is that region of the upper atmosphere where charged par- 
ticles (electrons and ions) of thermal energy are present, which are the result 
of ionization of the neutral atmospheric constituents by electromagnetic and 
corpuscular radiation. The lower boundary of the ionosphere, which is by no 
means sharp, coincides with the region where the most penetrating radia- 
tion (generally, cosmic rays) produce free electron and ion pairs in numbers 
sufficient to affect the propagation of radio waves (the D-layer). The upper 
boundary of the ionosphere is directly or indirectly the result of the interac- 
tion of the solar wind flow outward from the Sun at supersonic speeds with 
the planet. For weakly or essentially non-magnetic planetary bodies (Venus, 
Mars, Saturn’s satellite Titan), the interaction region between the solar wind 
and the ionospheric plasma represents the termination of the ionosphere on 
the sunward side; it is called the ionopause. 
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On the nightside the ionosphere can extend to greater distances in a tail- 
like formation, representing the solar wind shadow. In the tail the extent 
of the ionosphere is limited by the condition for ion escape. For magnetic 
planets with atmospheres (Earth, Jupiter, Saturn, Uranus, Neptune), the 
ionosphere terminates within the magnetosphere which comprises all charged 
particles of low (thermal) and high energies (radiation belts) . In this case the 
solar wind interacts with the intrinsic planetary magnetic field terminating 
at the magnetopause. The termination of the ionosphere is then the indirect 
result of the solar wind interaction; e.g., in the case of Earth, represented 
by the boundary between solar-wind induced convective motions inside the 
magnetosphere and the corotating ionospheric plasma called the plasmapause. 
(The region inside the plasmapause is also called the plasmasphere ; however, 
according to our definition this is simply part of the ionosphere.) 

1.1 Extreme Ultraviolet Radiation (EUV) and X-Rays 

The extreme ultraviolet (EUV) radiation nominally spans the wavelength 
range from 1000 to 100 A, although the edges are often somewhat indistinctly 
defined and extend shortward into the soft X-ray or longward into the far 
ultraviolet (FUV). Like X-ray emission, the production of EUV photons is 
primarily associated with the existence of hot gas in the Universe. 

Sources of EUV radiation can be divided into two main categories [1]. The 
first is where the emission arises from recombination of ions and electrons in 
a hot, optically thin plasma, which gives rise to emission line spectra; and 
the second is where objects are seen by thermal emission from an optically 
thick medium, resulting in a strong continuum spectra which contain features 
arising from transitions between different energy levels or ionization stages 
of several elements. 

Single stars and binary systems containing active coronae, hot O and 
B stars with strong stellar winds, supernova remnants and galaxy clusters 
belong to the former category. Hot white dwarf stars, central stars or plane- 
tary nebulae and neutron stars are possible continuum sources. Binary stars, 
where material is being transferred from a normal main sequence star, which 
is usually an M dwarf star onto a white dwarf, may well contain regions of 
both optically thin and optically thick plasma [1]. 

The emission of EUV radiation from astrophysical objects can also be used 
as probes of the interstellar medium (ISM). The absorbing effect of the ISM 
can modify the flux received at the Earth allowing, if the radiating source is 
well understood, its structure and density to be studied. Further, X-ray and 
EUV observations of various solar-like low mass stars (F, G, K, M) shown in 
Fig. 1.2 can be used for studying stellar-planetary relations, the ionization 
and heating of upper atmospheres of short periodic giant exoplanets and 
more interestingly “hypothetical” terrestrial exoplanets, which are expected 
to be discovered in the near future. The main characteristics of stars are 
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Table 1.1. Characteristics of stellar spectral types. 



Star type 


Main characteristic 


Color 


O 


Lines of highly ionized atoms are present, 
and lines of He + are visible in the spectrum. 


“Blue” 


B 


Neutral He lines are strong. Lines of Mg + 
and Si + are present. Hydrogen lines are 
conspicuous. 


“Blue” 


A 


Hydrogen lines are very strong. Lines of once- 
ionized and neutral metals (e.g., Fe) are also 
present 


“White” 


F 


Hydrogen lines are weaker than in class A, but 
still clearly visible. The ion Ca + has 
conspicuous lines. 


“White” 


G 


The spectrum looks like the Sun’s spectrum. 
The lines due to Ca + are very strong. 


“Yellow” 


K 


Many strong lines due to neutral atoms. 
Some molecular bands are present. 


“Orange” 


M 


Many strong lines due to neutral atoms are 
present in the spectra of stars of class M. 
Molecules such as TiO have strong bands in 
the spectra. 


“Red” 



listed in Table 1.1. Until the early 1970s the conventional view was that 
EUV astrophysics was not a practical proposition, because most elements 
have outer electron binding energies in the range of 10 to 100 eV; photons 
in the corresponding energy range will be strongly absorbed in any photon- 
atom interaction when the ionization potential IP < f)c/A, where A is the 
photon wavelength, f) the Planck’s constant and c the speed of light. When 
10 < IP < 100 eV, A lies in the range 1000 to 100 A, i.e., within the EUV 
band [1]. As a result, a planetary atmosphere is opaque to EUV radiation as 
in the case of Earth due to photoabsorption by N 2 , O 2 and O, so that the 
1/e absorption depth (unit optical depth) at 100 A is at an altitude of about 
120 km. 

The source of the EUV (1750 -70 A) and X-ray radiation (170-1 A) of 
a solar-like G star such as the Sun is produced in the various layers of the 
stars’ atmosphere. Table 1.2 lists the appropriate wavelengths at their point 
of origin [2, 3]. 

The absorption of the XUV radiation in a planetary atmosphere leads to 
photoionization and photodissociation of atmospheric constituents. Table 1.3 
lists the penetration depth of the solar radiation and particles as well as their 
energy fluxes. 

The relevant wavelengths for the heating of upper atmospheres are the 
ionizing wavelengths less than 1000 A, which contain only a small fraction 
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Fig. 1.2. Hertzsprung-Russell diagram: stellar bolometric luminosity in solar units 
as a function of spectral type of various stars. One can see that low mass K and M 
stars are more numerous than solar-like G-type stars. 

Table 1.2. XUV wavelengths and their origin in the solar atmosphere. 



A [A] 


Source of radiation 


>1800 


Photosphere, 


1200-2000 


Transit ion photosphere-chromosphere , 


900-1800 


Chromosphere, 


100-1000, Lya (1216) 


Transition to corona 


10-200 


Quiet corona, 


5-100 


Coronal active region 


1-50 


Thermal radiation . 




> solar flares 


0.01-10 


Non- thermal burst ) 



10 A= 1 nm 



[4] of the present solar spectral power as shown in Fig. 1.3 and Table 1.3. 
Figure 1.3 shows the altitude of penetrating solar radiation on Earth as a 
function of wavelength. 

The molecular constituents H 2 O, CO 2 and O 2 can be dissociated by rel- 
atively long wavelength radiation in the Schumann-Runge continuum range 
of 1750-1300 A, according to 

0 2 + l)i/(A < 1750 A) -> 0( 3 P) + 0( 2 D) 



C0 2 + l)i/(A < 1670 A) CO + OCD) . 
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Fig. 1.3. Penetration depth of solar radiation as a function of wavelengths in the 
Earth’s atmosphere. D, E and Fi refer to peaks of ionospheric Chapman layers. 
Ionization thresholds are also indicated. 



Table 1.3. Solar radiation and its penetration in the Earth’s atmosphere. 



A [A] 


Penetration depth 


Energy flux [W m 2 ] or [10 3 erg cm 2 sec x ] 


> 3000 visible + IR 


Troposphere 


1360 (outside the atmosphere) 


2000-3000 


Stratosphere 


14 


UV 


Mesosphere 




1500-2000 


Mesosphere 




EUV 


Thermosphere 


0.25 




Ionosphere 




1000-1500 


Thermosphere 


CO 

1 

O 

T— 1 

X 

00 


EUV 


Ionosphere 




100-1000 


Thermosphere 


2x 10 -3 


XUV 


Ionosphere 




10-100 


Thermosphere 


5xl0 -5 


X-rays 


Ionosphere 




1 


Mesosphere 


10“ 4 


X-rays 


Ionosphere 




Particles 


Solar protons 


Stratosphere 


10“ 3 


(MeV) 


Mesosphere 




Solar wind 


Magnetosphere 


10“ 4 


(keV) 
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The wavelength range with an energy flux of > 2 erg cm -2 sec -1 represents 
the predominant XUV heat source for the terrestrial thermosphere. Dissoci- 
ation can also take place below the true dissociation limit due to excitation 
of molecules into a state which dissociates. This predissociation occurs for 
O 2 , CO 2 and H 2 , but not significantly for N 2 , accounting for its apparent 
stability against photodissociation. 

The dissociation of O 2 is the first step in the old Chapman process of 
ozone (O 3 ) formation by a three-body collision 

o + o 2 + m— » 0 3 + M 

where M is a molecule ensuring momentum balance. The photochemical loss 
of O 3 occurs via dissociation and combination with O 

O + O 3 O 2 + O 2 . 

More recently catalytic loss processes for O 3 have been identified according 
to the following scheme 



o 3 + x xo + 0 2 

XO + O X + 0 2 

O 3 + O — > O 2 + O 2 

where the catalyst X are trace gases like NO, Cl, Br but also H and OH. 
O 3 absorbs in UV in the range 3000-2000 A and is therefore important for 
damage of DNA, absorbing most strongly at 2600 A. 

For H 2 , radiation A < 850 A leads to either dissociation or ionization. 
Atomic hydrogen possesses a strong continuum absorption cross-section be- 
low 912 A and strong Lyman-series lines, principally the intense Lyman-o; 
line at 1216 A. 

The ionization thresholds (ionization potentials, IP) for common con- 
stituents of planetary atmospheres are listed in Table 1.4 (Note: A [A] 
= 12,395/IP [eV]). Accordingly, XUV radiation of wavelengths below the 
Lyman-ce line (1216 A) is primarily responsible for the formation of plane- 
tary ionospheres and thermospheric heating. However, trace constituents, of 
meteoric origin, can be ionized by even longer wavelengths. Typical exam- 
ples are: Na (2410 A/5.14 eV); Ca (2028 A/6.11 eV); Mg (1622 A/7.64 eV); Fe 
(1575 A/7.87 eV); and Si (1520 A/8.15 eV). 

Rocket and satellite observations have led to a detailed identification of 
the solar emission line spectrum; relatively accurate values for the intensities 
of these radiations have become available [5, 6 , 7]. Solar photon fluxes at 
Earth for moderate activity are summarized in Table 1.5 [ 6 ]. The solar photon 
flux <£00 is used for the ionization process, while the XUV energy flux I 0 0 = 
(fyz/^oo is applicable to thermospheric heating. 

The solar photon fluxes vary both over a long period (11 year solar cycle) 
as well as over short periods (27 days) during disturbances. Although the 
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Table 1.4. Ionization thresholds/ionization potentials of atmospheric constituents. 





Constituent 


A [A] 


IP [eV] 




NO 


1340 


9.25 




nh 3 


1221 


10.15 




C 


1100 


11.3 




o 2 


1026 


12.1 




H 2 0 


985 


12.60 




ch 4 


954 


13.00 




H 


912 


13.59 




0 


911 


13.61 




C0 2 


899 


13.79 




CO 


885 


14.0 




h 2 


804 


15.41 




n 2 


796 


15.58 




A 


787 


15.75 




Ne 


575 


21.56 




He 


504 


24.58 


Table 1.5. Solar XUV photon and energy fluxes. 


A [A] 


#oo [10 9 ph cm 2 sec 


1 ] loo [erg cm 2 sec 1 


1215.7 (Ly a) 


300 




5 


1027-911 


11.61 




0.23 


(1025.7, Ly 0) (3.5) 




(0.067) 


(977, C III) 


(4.4) 




(0.090) 


911-800 


8.3 




0.20 


800-630 


2.4 




0.064 


630-460 


4.7 




0.17 


(584.3, He I) 


(0.9) 




(0.03) 


460-370 


0.63 




0.03 


370-270 


10.3 




0.65 


(303.8, He II) 


(5.4) 




(0.35) 


270-205 


4.5 




0.36 


205-153 


4.6 




0.49 


153-100 


0.4 




0.06 


120-80 


0.3 




0.066 


80-40 


0.33 




0.108 
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Table 1.6. Variations in the X-ray energy flux 7x with solar activity. 



Condition of the Sun 


7x [erg cm 2 sec 1 


1 

2 A 


4 A 


6 A 




Completely quiet 


10 -8 


10- 7 


icr 6 


Solar minimum 




Quiet 


10" 7 


10" 6 


10“ 5 




Lightly disturbed 


10“ 6 


10" 5 


10 -4 


Solar maximum 




Disturbed 


10“ 5 


10“ 4 


10“ 3 




Special events 


10“ 4 


10“ 3 


10“ 2 


Flares 




Class 3 flares 


icr 3 


10“ 2 


lO -1 



exact amplitude of these variations is not yet fully established, it appears to 
be for the EUV range of the order of 2. The 10.7 cm solar radio flux F10.7 is 
usually considered to be an excellent indicator of solar activity; according to 
limited data the integrated flux in the range 1310-270 A changes by a factor 
of about 1.5 for a change in F10.7, by a factor of 2 [5]. 

The solar cycle variation for the range 100-10 A amounts to a factor of 7; 
the energy flux in this range is ~0.9ergcm -2 sec -1 at solar maximum. The 
X-ray energy flux /x at 8 A in the absence of flares varies over the solar cycle 
by a factor of 300; the variation in the range 8-2 A from a completely quiet 
Sun to a class 3 flare can amount to five orders of magnitude (10 5 ), as shown 
in Table 1.6 [3]. 

Production of Ionization by Photons 

In the photoionization process the production of an ion pair (ion 4- electron) 
occurs as the result of the action of a photon whose energy \)v is equal to or 
greater than the ionization potential of the constituent X, according to 

X + |)i/(> IP) -> X + + e . 

For the general discussion of the production of electrons and ions by pho- 
toionization of neutral atmospheric constituents we will make some simpli- 
fying assumptions. The more general case can easily be obtained from this 
starting point by summation over all wavelengths. 

The following derivations are based on the assumptions of 1) monochro- 
matic radiation, i.e., a small wavelength interval over which the absorption 
and ionization cross-sections can be taken as constant; 2) an isothermal at- 
mosphere where the atmospheric scale height is H— const. ( H — k^T/mg: 
cf. Sect. II. 1) and a horizontally stratified atmosphere (i.e., neglecting the 
curvature of the planetary surface). The ionizing radiation is assumed to be 
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incident at an angle \ with the vertical (the solar zenith angle), having a 
photon flux outside the atmosphere $oo. The ion production rate depends 
on the number density of the ionizable constituent n j, its ionization cross- 
section G\ for the particular wavelength interval and the local photon flux 
according to 



q = . 



(i.i) 



Along a path length d s, nj<r a ds photons are absorbed, where cr a is the ab- 
sorption cross-section which is usually larger than a\. Hence we can write 

d<Z> 

— = rijcr a^Poc sec% 
dz J 

where ds = dz sec x, and obtain by integration 



& = ^oo exp 



-<r a secx 



POO 

J n A z ) 



dz 



( 1 . 2 ) 



i.e., the ionizing photon flux decreases from its value outside the atmosphere 
as the result of absorption. The ion production rate can now be expressed for 



n j (z)=n j0 ex p (—g) 



as 



q = crjnjo^oo exp 



* r , ' 

— - cr a sec x y naz 



(1.3) 



or in terms of the optical depth r = cr a sec x J^° n j dz = a a n(z)H sec x 



q = airij(z)4> 00 e r . (1.4) 

The ion production rate has a maximum where dq/dz = 0, which also cor- 
responds to the condition r = 1, for overhead Sun. Accordingly, the altitude 
where the maximum of ion production occurs can be found from this condi- 
tion. Since 



r = a a rijoH sec x * exp 



H 



= 1 



we have exp (z*/H) = a a noH • sec% where it should be noted that z* is the 
relative height above the reference level z = 0, represented by an absolute 
altitude ho, usually close to or above the homopause. 

Hence the absolute altitude of the ion production maximum, h*, is given 
by 



h* = H \n(a a noH sec x) + ho = + H In sec x • 



(1.5) 
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The maximum value of the ion production rate is given by 

$oo 

Qm = — cos X = q 0 cos x 

o a eH 



( 1 . 6 ) 



where the ratio cr* /cr a = r\i is called the ionization efficiency. The condition of 
r = 1 for overhead Sun (x = 0) is often used to indicate the penetration depth 
of ionizing radiation, i.e., where it decreases by 1/e or where for H =const., 
the ionization rate would reach an absolute maximum 

Vi^oo 

qo = ~^r 

The altitude Hq is also called the altitude of unit optical depth , i.e., where 
a a n(h* 0 )H = 1. It should be noted that the height of maximum ion production 
depends on the optical depth and solar zenith angle (secx), but not on the 
photon flux, whereas the maximum value of the ion production rate depends 
on the photon flux outside the atmosphere, the scale height and cosx- 

The ion production rate as a function of altitude can be written using the 
above definitions 



Q = 



q m exp 




h-h* 

H 



— sec x exp 



(h-h*) 



H 



} 



(1.7) 



This is the well-known Chapman ion pair production function [8] which is 
illustrated in Fig. 1.4 in normalized form 



q = qo cosxexpjl — z — secxexp(— z)} (1.8) 

where z = (h — h$)/H. The optical depth and thus the penetration of ionizing 
radiation depends on the absorption cross-section. In the XUV range, most 
major constituents of planetary atmospheres have absorption cross-sections 
of order a a ~ 10 -17 -10 -18 cm 2 , whereas at X-ray wavelengths (see Table 1.7) 
the cross-sections are much smaller [9, 10]. Between 50 and 20 A the K ab- 
sorption limits occur, which correspond to an abrupt change in cr a by about 
one order of magnitude. The (total) absorption cross-sections are the upper 



Table 1.7. 



A [A] 


<7 a [cm 2 ] 


1 


~ 10“ 22 


3 


~ io- 21 


5 


~ 10“ 20 


15 


~ 10- 19 


50 


~ 10“ 19 


100 


00 

1 

0 

rH 

1 
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Fig. 1.4. Ion-pair production function, as a function of solar zenith angle x> nor- 
malized to the absolute ionization maximum q for overhead Sun [x = 0] occurring 
at an altitude h q. The height parameter corresponds to z = (h — ho)/ H. 



limit for the ionization cross-sections; for atomic species cr a = cq, i.e., the ion- 
ization efficiency is r][ = 1, whereas for molecular species o\ < cr a , i.e., 77 * < 1. 

Figs. 1.5-1. 8 show the absorption cross-sections for various important 
constituents of planetary atmospheres [9]. 

The ion production rate derived above under a number of simplifying 
assumptions can be extended to be more generally applicable. For an at- 
mosphere with a constant scale height gradient dH/dz — (3, the ionization 
maximum occurs where r = 1 + f3 and the ion production rate is given by 

q = <?m(l + /?){! — C — e _c sec \} 

where ( is a height variable defined by H = Hoe ^ . 

Near sunrise and sunset, i.e., x > 75°, the sec x term in the ion production 
function, which is a consequence of the flat planet approximation, has to be 
replaced by the so-called Chapman function Ch(x) [11]. 



Ch(x) = 



J s °° n(s) ds 

JT n ( h ) dh 



This function is defined as the ratio of the total content n of the atmosphere 
in the line of sight (5) to the Sun to the vertical content ( h ). 





16 I Radiation and Particle Environments from Stellar and Cosmic Sources 




Fig. 1.5. Absorption cross-section of N 2 from experimental data (left); O 2 , (right) 
(after E. W. McDaniel [9]). 





Fig. 1.6. Absorption cross-section of O (left) and N (right) from theoretical calcu- 
lations (after E. W. McDaniel [9]). 




Fig. 1.7. Absorption cross-section of H 2 based on experimental data (after E. W. 
McDaniel [9]). 
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X (A) 



Fig. 1.8. Absorption cross-section of CO 2 from experimental data (after E. W. 
McDaniel [9]). 

A number of analytical approximations to the Chapman function have 
been developed, including the extension to a constant scale height gradient (3 
(cf. Sect. II. 1). For an isothermal atmosphere a useful approximation of the 
Chapman function is given in terms of the tabulated error function 




where x = R/H — (Rq + h)/H with H the scale height and Rq the planetary 
radius. The Chapman function for different values of x and its comparison 
with sec% is shown in Fig. 1.9. 

For x = 90°, Ch(x, 7r/2) = (x7r/2)^, representing the ratio of the total 
content in the line of sight to the vertical content of an atmospheric con- 
stituent, a quantity of interest in occult at ion experiments. At high enough 
altitudes where the optical depth is small the exponential factor in (1.9) can 
be neglected and the ion production function in this “low attenuation region” 
assumes the simple form 



q = J(a u \)nj(z) (1.10) 

where J = ( 01 ^ 00 ) is the photoionization rate coefficient (or frequency of 
ionization ^ on , [sec -1 ]). The ionization coefficient is based on the averaged 
ionization cross-sections and photon fluxes over all pertinent wavelengths. For 
the light atmospheric constituents H and He, which will make up the main 
constituents of planetary exospheres, the ionization coefficients are listed for 
Venus, Earth and Mars in Table 1.8. Since many wavelengths will contribute 
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Fig. 1.9. Chapman function Ch(x, x) for parametric values of x — R/H , compared 
with sec x- 



Table 1.8. Photoionization rate coefficients. 



J [sec x ] 


Venus 


Earth 


Mars 


Atomic hydrogen (H) 
Helium (He) 


2 x 10" 7 
10“ 7 


1.5 x 10“ 7 
6 x 10“ 8 


4 x 10 -8 
2.5 x 10“ 8 



to the ionization of a particular constituent rij and will also be absorbed by 
other atmospheric constituents, the ion production rate in this generalized 
form can be expressed by 

nj<Tij$x uv(A)exp 

In the photoionization process, photoelectrons are produced whose energy, 
E pe = fyv — IP, may be large enough, especially in the X-ray range, to lead 
to additional (“secondary”) ionizations. In the X-ray range, the ion produc- 
tion is also computed by using the X-ray energy flux lx divided by 35 eV, 
the average energy expended in the formation of an ion pair for air or CO 2 . 
It should be noted that the calculation of secondary ionization depends on 




£ 

l 






rij sec x dz 



(1.11) 
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the initial photoelectron distribution and their collision processes with atmo- 
spheric particles. 

1.2 Corpuscular Radiation 

Galactic cosmic rays, energetic solar particles, solar wind protons and lower 
energy electrons and protons populating a planetary magnetosphere are re- 
sponsible for the formation of electron-ion pairs at the lower level of planetary 
ionospheres [12]. Because of their high energy, cosmic rays penetrate deepest 
into the planetary ionosphere. 

Galactic cosmic rays with energies E > 10 9 eV are responsible for the 
formation of the lowermost ionospheric layer (sometimes called the cosmic 
ray or C-layer) actually representing the lower part of the D-layer. 

Relativistic solar cosmic rays (protons) produced during highly disturbed 
solar conditions, which are a rather rare occurrence (once in several years) 
having energies E > 10 BeV, cause a substantial enhancement in the ioniza- 
tion of the D-layer [13]. 

Subrelativistic solar cosmic ray protons with energies 10 MeV < E < 
100 MeV produce strong ionization effects at polar caps which are accessible 
on a magnetic planet, leading to polar cap absorption (PC A) events in the 
terrestrial ionosphere. These events are generally observed at periods of max- 
imum solar activity, but are almost non-existent during solar minimum [14]. 
The general formula for ion-pair production by charged particles (corpuscular 
radiation) can be written [15] 

q(h) = ^Un^~ i{E)iEAa (112) 

where W = 35 eV is the average energy required for the formation of an 
ion pair in air or CO 2 , dE/dx = secxdE/dh is the energy loss in an in- 
elastic collision process (ionization loss), j(E) = KE ~ 7 represents the dif- 
ferential energy spectrum of the ionizing particles and i? is the solid angle. 
j(E) =dJ/dE , where J(> Eo) = f£?j(E)dE is the integral (spectrum) flux 
for energies grater than Eo. 

In the presence of a planetary magnetic field, the penetration of charged 
particles is inhibited; for a dipolar field only the most energetic particles can 
penetrate to low latitudes, so that the particle ionization effects are normally 
restricted to high magnetic (polar) latitudes. The ability of a charged particle 
with momentum (mv) to penetrate a magnetic field of strength B can be 
expressed by the rigidity P which is defined by 

TYIVC 

P — —— [Volt] = Btb [Gauss/cm] 

Ze 

where Z is the charge number, e is the electronic charge, and vb is the Larmor 
or gyro radius (cf. Chap VI) of a charged particle of momentum mv moving 
in a plane perpendicular to the magnetic field vector B. 
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The ion production rate resulting from galactic cosmic rays can be written 
in simplified form [13] 



q(<p) = 



go (<p)n 

no 



[cm 3 sec x ] 



(1.13) 



where qo(ip) is the ionization rate at magnetic latitude (p for an atmospheric 
level with the corresponding number density no = 2.6 x 10 19 molcm -3 . The 
ionization rate at geomagnetic latitudes ip > 70° corresponds to that which 
would prevail if the galactic cosmic rays were unimpeded by a planetary 
magnetic field; the ionization rate at the geomagnetic equator is reduced 
by an order of magnitude relative to that at the pole. The polar ionization 
rate at solar maximum has been found to be go = 200 cm -3 sec -1 atm -1 , 
whereas that at solar minimum is larger by about a factor of 4, while that 
at middle and low latitudes varies by about a factor of 2, ranging from 150 
to 300 cm -3 sec -1 at (p = 50° [13]. This inverse relationship of cosmic ray 
intensity with solar activity, also called the solar modulation of cosmic rays , 
is due to the interplanetary magnetic field carried by the solar wind [14]. 
The penetration depth h of charged particles in a planetary atmosphere is 
expressed by the range 91 in units of gem -2 , according to 

poo 

91 = sec x / p(h)dh (1.14) 

Jh 



where x is the zenith angle of the incident particles, and p is the mass density 
(p = nm). 

The range 91 is related to the ionization loss or stopping power according 
to 



91 



-j 1 . 

J E(h) 



1 cL E 

[ p d h 



i -i 



d E. 



(1.15) 



Range-energy relations are usually based on experimental data or theoretical 
expressions based on the Bethe formula [16]. Empirical formulas exist for air, 
and since the pertinent physical parameters for CO 2 are very close to those 
for air, they may be applied to the atmospheres of the terrestrial planets. 
The range of protons with energies E < 500 MeV is approximately given by 
[17] 



Dip [gem 2 ] = 



E 1 - 78 

420 



and for electrons of energy E < 200 keV by 



91e[gcm 2 ] 



£ 1.96 

0.75 



where E is the particle energy in MeV. 
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Since charged particle trajectories are not rectilinear due to changes in 
direction in the collision process, as well as due to a planetary magnetic 
field, accurate analysis of energy deposition has to take these effects into 
account. Such calculations have been performed for low energy electrons using 
a statistical Monte Carlo method [18]. 

The approximate penetration depths at vertical incidence of energetic 
electrons and protons in the atmospheres of Earth, Venus and Mars, based 
on range-energy relations, are shown in Fig. 1.10. The most penetrating 
charged particles are galactic cosmic rays and energetic solar protons, which 
produce ionization at altitudes below 100 km while lower energy particles, 
such as auroral electrons (1-10 keV), deposit their energy at altitudes 100km 
in the terrestrial atmosphere. Charged particles penetrate to lower altitudes 
on Mars, than on Earth, but not as deep into the atmosphere of Venus as on 
Earth. This is a consequence of the fact that in the lower thermospheres of 
these planets the scale heights are not greatly different, whereas the densities 
are highest for Venus and lowest for Mars, with Earth lying in between. On 
Saturn’s satellite Titan an ionospheric layer produced due to cosmic rays is 
expected at an altitude of about 90 km. 

For energetic solar protons and galactic cosmic ray protons, the use of 
an exponential rigidity spectrum of the form J(> P ) = J 0 exp(— P/P 0 ) with 
Po < 100 MeV for polar or non-magnetic conditions and P 0 > 150 MeV 
for geomagnetic latitudes p < 65°, has been found to give a good rep- 
resentation of their ionization rate [13]. Fig. 1.11 shows the ion-pair pro- 
duction functions over the polar cap for galactic cosmic rays during solar 
maximum and minimum and a solar proton event with a rigidity spectrum 
J > 100 MV = 0.05cm _2 sec _1 ster -1 and Po = 50 MeV [13]. (Some solar 
proton events result in a peak ionization rate 100 times larger than this.) 



ELECTRON ENERGY {k*V) 




1. 10 Approximate altitudes of penetration of energetic particles in the atmospheres 
of Earth (®), Mars (cf) and Venus ($). 
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Fig. 1.11. Ionization rate due to galactic cosmic rays and high energetic solar 
protons in the terrestrial atmosphere over the polar cap (after W. Webber [13]). 



The solar proton and galactic cosmic ray ionization rate for Mars and Venus 
[13, 18, 19] can be scaled from the terrestrial case according to 






p(h m , v ) 
p(h E ) 



q{h e ) 



(1.16) 



where q(h e) is the ionization rate at the terrestrial poles, He and hu,v 
are the atmospheric penetration altitudes resulting from the appropriate 
range 91(gcm -2 ) for the same cutoff energy for solar protons, and pQie) 
and p(hy^M) are the respective mass densities. For low energy electrons and 
protons, such as auroral particles, the ion production rate due to particles in 
the energy range E m - m to 2£ max can be expressed for Earth by 

^^C“ L ^ E ’ h)iE < li7 > 

u -C'min 

where p(h) is the atmospheric density in units of gem -3 at altitude h. The 
ionization loss for electrons can be approximated by [17] 

L e (E) = — ^ = 0.38£T°- 96 
v ’ d« 

and for protons by 

L p ( E ) = -^ = 236E-°- 78 



where 91 is the range in g cm 2 and E is in MeV. 
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The number of particles of energy E arriving at altitude h , i.e., the dif- 
ferential energy spectrum j(E , h) in cm -2 sec -1 ster -1 keV -1 , has to be ob- 
tained by integrating over both pitch (a) and solid (i?) angle if a magnetic 
field is present, according to 



j{E,h)=f j j(E) exp 
Jo Jo 



X (E) cos a 



sin a da d i? 



(1.18) 



where x is the range at height h and A (E) — 3.15 x 10~ r E 22 is the attenuation 
mean free path. 

Another possible ionization source for planetary ionospheres is the solar 
wind. Although most of the solar wind flows around the planetary obstacle 
(cf. Chap. VI), a small fraction of the solar wind protons may leak through 
the magnetopause or ionopause, especially in the tail region. 

Ionization by the stellar wind is due to two processes: 1) direct impact ion- 
ization of atmospheric constituents by stellar wind protons having an energy 
E ~ 1 keV, and 2) ionization by hot hydrogen generated by charge exchange 
(cf. Chap. IV) between stellar wind protons and atmospheric constituents. 
The ionization rate due to the stellar wind can be expressed by 



q = a' i n j ^' sw + 

where the g\ are the ionization cross-sections (~ 10“ 11 cm 2 ), #' w is the stellar 
wind proton leakage flux, i.e., a small fraction (< 10%) of the stellar wind 
flux outside the planetary magnetopause or ionopause ^ sw , and is the 
effective ionizing flux of hot hydrogen. The latter is governed by a continuity 
equation of the form 



d^H* 1 d<P' sw d<2>H* dE 

ds 2 ds dE ds 

where the factor | denotes the downward part of an isotropic flux of 
hot hydrogen generated by charge exchange between stellar wind protons 
and the neutral constituent n^, according to d<£' sw /ds = — < 7 c n^ sw , with 
<r c « 10~ 15 cm 2 ; the second term on the right-hand side represents the ther- 
malization of the initially formed hot hydrogen beam along its path by scat- 
tering and ionization processes (cf. Chap. III). 



1.3 Evolution of the Solar/Stellar XUV Flux with Time 

Because the type of the star plays a major role in all atmospheric processes, 
the evolution of planetary atmospheres must be understood within the con- 
text of the evolving stellar energy and particle fluxes [4]. Because thermal 
escape of atmospheric constituents in planetary atmospheres depends on the 
stellar XUV radiation, which affects the temperature of the upper atmosphere 
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(cf. Chap. II) and photoionization processes, it is important to study the time 
evolution of stellar radiation and particle fluxes. 

Observations of enriched heavy isotopes in planetary atmospheres [20, 21] 
and radiative stellar fluxes, stellar magnetic fields, stellar winds of solar-type 
stars with different ages [22, 23, 24, 25] and lunar and meteorite fossil records 
[26] indicate that our early Sun underwent a highly active phase after its 
arrival on the zero-age main-sequence (ZAMS). 

Solar-Like G Stars 

The observational evidence of an active XUV phase of the young Sun is also 
a representative for other solar-like G stars and lasted about 0.5 to 1.0 Gyr. 
It included continuous flare events, where the radiation environment was 
several hundred times more intense than today. The high radiation levels of 
young solar-like stars is triggered by strong magnetic activity. The magnetic 
activity is expected to have greatly decreased with time [24, 27, 28] as the 
solar rotation slowed down through angular momentum loss. Observational 
evidence and theoretical models suggest that young solar-like stars rotate 
about 10 times faster than the Sun today and have significantly enhanced 
magnetically generated coronal and chromospheric activity [24, 29]. 

Solar proxies are studied inside the so-called Sun in Time program [24], 
which was established to investigate the magnetic evolution of the Sun by 
using a homogeneous sample of single nearby GO-V main sequence stars, 
which have known rotation periods and well-determined physical properties, 
including temperatures, luminosities, metal abundances and ages. As can be 
seen in Table 1.9, the sample of studied solar proxies contains stars that cover 
most of the Sun’s main sequence lifetime from 130 Myr to 8.5 Gyr. 

A large amount of multiwavelength (X-ray, EUV, FUV, UV, optical) data 
have already been collected. The observations of solar proxies, obtained with 
the ASCA, ROSAT, EUVE, FUSE and IUE satellites, cover a range between 
3300 and 1 A, except for a gap between 920 and 360 A, which is a region 
of very strong interstellar medium absorption. Details of the data sets and 
the flux calibration procedure employed are provided in [24]. Full spectral 
irradiance tables have already been completed for five of the stars in the Sun 
in Time sample (EK Dra (130 Myr), 7T 1 UMa (300 Myr), k, 1 Cet (750 Myr), 
/ 3 Com (1.6 Gyr), (3 Hyi (6.7 Gyr)) and show an excellent correlation between 
the emitted energy flux and stellar age. 

The coronal XUV emissions of the young main-sequence Sun were about 
100 to 1000 times stronger than those of the present Sun. Similarly, the 
transition region and chromospheric FUV-UV emissions of the young Sun 
are expected to be 10 to 100 and 5 to 10 times stronger, respectively, 
than present and the flux variation over age is therefore a steep wavelength 
function. Figure 1.12 shows the time evolution of the spectral range with 
1 A < A < 1000 A, which includes X-rays and EUV and the Lyman-a line at 
1215.6 A [30] at a distance of 1 AU. The resulting relative XUV fluxes shown 
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Fig. 1.12. Time evolution of the Ixuv energy flux for solar-like G stars (solid line: 
A = 1000-1 A; dashed line: Lyman-a = 1215.6 A). 



in Table 1.9 yield an excellent correlation between the emitted flux and stel- 
lar age. In the 1000-1 A interval, the fluxes follow a power-law relationship 
[30] 



/xuv(t ) = 6.16 x (t[Gyr])~ 1-19 . (1.19) 



Table 1.9. Solar-like G-type stars studied within the Sun in Time program. The 
parameters of the solar proxies shown below are: luminosity, distance, age, the 
XUV energy flux relative to the present value at 1 AU in the 1000—1 A interval 
and Lyman-o [24]. 



Star 


Lum [L©] 


Dist. [pc] 


Age [Gyr] 


Prot [d] 


UcuvO) 

Cxuv 


J Lq(«) 

J L a 


EK Dra 


0.94 


34 


0.13 


2.75 


69.8 


14.6 


7r 1 UMa 


0.98 


14.3 


0.3 


4.68 


25.8 


7.8 


Xl Ori 


1.07 


8.7 


0.3 


5.08 


25.8 


7.8 


9 Cet 


0.98 


20.4 


0.65 


7.6 


10.3 


4.4 


k1 Cet 


0.84 


9.2 


0.75 


9.2 


8.7 


3.9 


0 Com 


1.36 


9.2 


1.6 


12.4 


3.5 


2.2 


15 Sge 


1.25 


17.7 


1.9 


13.5 


2.9 


1.9 


Sun 


1.00 


- 


4.6 


25.4 


1 


1 


18 Sco 


1.05 


14.0 


4.9 


23.0 


0.9 


0.96 


0 Hyi 


3.58 


7.47 


6.7 


28.0 


0.6 


0.76 


16 Cyg A 


1.62 


21.6 


8.5 


35.0 


0.5 


0.63 





26 



I Radiation and Particle Environments from Stellar and Cosmic Sources 



At longer wavelengths, the Lyman-a emission feature can contribute to a 
significant fraction of the XUV flux. High-resolution Hubble Space Telescope 
(HST) spectroscopic observations were used to estimate the net stellar flux. 
These measurements, together with the observed solar Lyman-a, define the 
following power-law relationship with high correlation 

= 3.17 x (f[Gyr]) -0 ' 75 . (1.20) 

In both power laws, the XUV and Lyman-a expressions are valid for ages 
between 0.1 and 7 Gyr, /xuv and /l q are the present integrated fluxes at 1 
AU and /xuv(t) and /L a (t) are the integrated fluxes as a function of time. 
One finds fluxes of « 6 x /xuv and « 3 x Ii Ja about 3.5 Gyr ago, and 
« 100 x /xuv and « 20 x /l q about 100 Myr after a solar-like G-type star 
arrived on the ZAMS. 

Low Mass K and M Stars 

Similar studies of the properties of X-ray and extreme ultraviolet fluxes and 
their typical evolution for low mass K and M stars are important, since 
“hypothetical” terrestrial exoplanets may also orbit in low mass star systems. 

For an initial estimate of the evolution of XUV irradiances one can use a 
proxy indicator, which is the ratio of the X-ray luminosity to the bolometric 
luminosity (log [Lx/Lboi])- This ratio is highest for the more active stars (i.e., 
fastest rotation period) and decrease monotonically with decreasing level of 
chromospheric activity [31, 32]. 

From the analysis of open cluster stars it is well established that all single 
late-type G, K and M stars spin down as they age, their activity decreases 
with time (see also P TOt in Table 1.9), and so does the ratio log(Lx/Lboi)* 
Another effect is that log(Lx/Lboi) does not increase up to values arbitrarily 
close to unity for very active stars. Instead, a saturation phenomenon occurs 
and no active star seems to reach higher levels than log(Lx/Lboi) ~ —3 
[33, 34, 35]. 

Although some authors have interpreted this phenomenon as 100% sur- 
face coverage with stellar active regions like flares, the saturation phe- 
nomenon is not yet completely understood [1]. Qualitatively, the evolution 
of log(Lx/Lboi) for a late-type star has a flat plateau from its arrival on the 
main sequence up to a certain age (the end of the saturation phase) and then 
decreases monotonically with age. 

To study the activity level of K and M stars one can compare their values 
of log(Lx/Lboi) and age with those of solar-like G-type stars. With the same 
underlying physical mechanism responsible for XUV emissions (magnetic dy- 
namo) and a supposedly similar spectral energy distribution, it is reasonable 
to assume that stars with similar values of log(Lx/Lboi) will also have similar 
log(Lxuv/Tboi) (he., 1 A< A <1000 A). One can use the XUV flux evolution 
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relationships for solar-like G stars and infer the XUV fluxes of K and M stars 
provided their X-ray fluxes are known. Astrophysical observations show that 
solar-like GO-V stars stay at saturated emission levels until ages of ~100 Myr 
and then their X-ray luminosity rapidly decreases following the power-law 
relationship discussed before as a function of age. 

Recent studies of K type stars show that they stay at saturated emission 
levels for a little longer time and then also decrease following a power-law 
relationship of a very similar slope. Interestingly, M0-M5 stars seem to have 
saturated emission levels up to 1 Gyr (and possibly longer) and then decrease 
in an analogous way to G and K stars [36]. 

These preliminary results indicate that early K stars and early M-type 
stars may have XUV irradiances that are about 3 to 4 times and about 
10 to 100 times higher, respectively, than solar-type stars of the same age. 
More accurate investigations are currently being carried out with an extended 
sample and a large variety of observational data [36] . 

EUV observations by the Einstein satellite on late G, K and M stars 
indicate that their main and median EUV luminosities of G, K stars broadly 
match the X-ray luminosity functions, while the M star EUV luminosities may 
be about a factor of 10 lower [1]. This probably arises from a high luminosity 
M dwarf tail in the X-ray function, but more research and observational data 
are needed. 

However, high radiation on saturation levels up to several 100 Myr and 
even Gyr of M stars may result in large temperatures in the upper atmo- 
spheres on terrestrial exoplanets, which may result in problems of their at- 
mospheric stability, evolution of their water inventories and planetary habit- 
ability (cf. Chap. II, Chap. VI, Chap. VII). 



1.4 Solar/Stellar Winds 

After the XUV radiation environment of stars, the evolution of the stel- 
lar wind mass flux as a function of age also has important implications for 
the evolution of planetary atmospheres. The average momentum flux of the 
present solar wind ions at 1 AU is about 1.5 x 10“ 8 g cm -2 s _1 , while the 
energy flux is about 2 x 10 11 eV cm -2 s -1 . To obtain the evolution of the 
stellar wind velocity v sw and stellar wind density n sw of solar-like stars, one 
can use the scaling for stellar mass loss provided by [25] and the scaling for 
the velocity developed by [26]. 

The stellar mass loss can be estimated by using recent indirect stellar wind 
observations [25]. Recent HST high-resolution spectroscopic observations of 
the H Lyman-a feature of several nearby solar-like G and K main-sequence 
stars have revealed neutral hydrogen absorption associated with the interac- 
tion between the stars’ fully ionized coronal winds with the partially ionized 
local interstellar medium (ISM). By modelling the absorption features formed 
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in the astrospheres of these stars one can provide empirically estimated coro- 
nal mass loss rates for G and K main-sequence stars [25]. 

It was found from a small sample of solar-like stars, where astrospheres 
can be observed, that stellar mass loss rates increase with stellar activity. 
The observations of younger solar-like stars indicate that these stars had 
much denser stellar winds than the Sun today. The correlation between mass 
loss and X-ray surface flux follows a power-law relationship, which indicates 
an average stellar wind density up to 1000 times higher than today during the 
first 100 Myr after the Sun reached the ZAMS. The mass loss as a function 
of the X-ray flux (j>x observed on G and K stars, when the coronal activity 
is at its maximum, can be written as [25] 



d M 

d t 



OC0x 



1.15 



( 1 . 21 ) 



Mass loss rates of cool main-sequence stars depend on their rotation periods 
^rot> which are in turn correlated with the star’s ages. According to [25] these 
relations can be expressed as 



4>x « P- t 2 ' 9 (1.22) 

and 

Protect- 0 ' 6 (1.23) 

respectively. From (1.21) and (1.22), one obtains a power-law relationship 
for the mass loss rate as a function of rotation period, dM/dt(P rot ). Here, a 
complication occurs. Equation (1.21) does look like a scaling for the mass loss, 
but what is really measured by [25] is the total ram pressure, i.e., the product 
of mass loss and stellar wind velocity. The mass loss given was obtained by 
assuming a constant stellar wind velocity v sw . For finding scalings for both 
n sw and v sw , one can write Mv sw rather than M in (1.21) [37]. Thus, 

^-v sw <x p- t 3 - 3 . (1.24) 

To obtain the time dependence, one has to insert a function for P TO t(t) into 
(1.24). In principle, (1.23) could be used, but for sake of consistency with the 
velocity scaling given below one can take the scaling derived by [26] 

P rot oc (l + f ) (1.25) 



with the time constant r = 2.56 x 10 7 yr. The time elapsed since the formation 
of the planetary system is denoted by t. For the present day solar system, 
t = 4.6 Gyr. Combining (1.24) and (1.25) it is possible to derive a power-law 
formula for the stellar mass loss as 
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d M 

“TT^sw CX 




(1.26) 



The stellar mass loss linearly depends on v svl and n sw : 



d M . 

— dTZsw^SW^p 



(1.27) 



where A is the stellar surface area and ra p the mass of the stellar wind protons. 
Thus, the average time-behavior of the stellar wind velocity of solar-like stars 
can be achieved by [26, 38] 



v sw — v* 




(1.28) 



Prom the mass loss formula (1.26) and with (1.27) one can determine the 
evolution of the average particle density [37] 



'fl'SW — 




(1.29) 



The proportionality constants are determined by the present-day conditions. 
With v sw — 400 km s _1 and n sw = 10 7 m -3 for t = 4.6 Gyr and at r = 1 AU 




Fig. 1.13. Time evolution of the observational based minimal and maximal stellar 
wind densities scaled to 1 AU (left scale: solid lines) and obtained from several 
nearby solar-like G and K main-sequence stars due to observations from H absorp- 
tion associated with the interaction between the stars’ fully ionized coronal winds 
with the partially ionized local ISM. On the right scale one can see the temporal 
evolution of the stellar wind velocity (right scale: dashed curve). More observations 
of early active stars with ages less than 700 Myr are needed to obtain a better 
picture of the mass- loss/activity relation (courtesy of T. Penz). 
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[39] one obtains v* = 3200 km s _1 , n* = 2.4 x 10 10 m -3 (density at 1 AU). 
The time constant is r = 2.56 x 10 7 yr [26]. For distances other than 1 AU, 
(1.29) is scaled with a 1/r 2 dependency. The minimum and maximum time 
variation of n sw (t) obtained from the observations of [25] and v S w(t) [26] at 
1 AU are shown in Fig. 1.13. 

One can see from Fig. 1.13 that the observational data for solar-like G 
and K stars suggest that more active stars should have a higher mass loss 
rate and solar wind number density. However, recent observations of the M 
dwarf star Proxima Cen and the RS CVn system A And (G8 IV + M V) 
are inconsistent with this relation and show lower mass loss rates. A recent 
observation of the £ Boo binary system shows also a lower mass loss rate, 
which is consistent with the mass loss rates previously found for Proxima 
Cen and A And (J. L. Linsky, private communication, 2004). 

The common feature of these three stars is that they are all very active 
in X-ray surface fluxes at about 10 6 erg cm -2 s _1 , which is about a factor 
30 larger than the todays Sun and corresponds to a time of about 700 Myr 
after solar-like young stars arrived at the ZAMS. These recent observations 
indicate the uncertainty of the early mass loss and solar wind estimations, 
because both stars in the £ Boo system are usual but very active G and K 
stars, while Proxima Cen and A And are different types than normal G and 
K stars. 

It seems possible, that there could be a high-activity cutoff to the mass- 
loss/activity relation in (1.21) [25]; therefore, more active young solar-like G 
and K stars with X-ray surface fluxes larger than 10 6 erg cm" 2 s _1 should 
be studied in the future. 
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II. 1 Barosphere: Distribution Laws 

Static Atmosphere 

The upper atmosphere, in fact the entire barosphere, can be characterized by 
its pressure and density distribution. The starting point for the derivation of 
this distribution is the hydrostatic equation 

dp=-gpdz (2.1) 

and the perfect gas law 



p = nk B T (2.2) 

where p is the pressure, g is the acceleration of gravity, p — nm is the mass 
density with n the number density, m ~ J2j n j m j/ n j * s the mean molecular 
mass, ks = 1.38 x 10 -16 erg/K is the Boltzmann constant , T is the absolute 
temperature and z is the height variable. 

Combining (2.1) and (2.2) and assuming first the simplest case of T, g 
and m =const. we obtain after integration 



(mg \ ( z 

P = PoeXP {-^T Z ) =PoeXp (-ff 

n = n 0 exp (--^) 
p = p 0 exp (“) 



(2.3) 



where H — k B T/mg ( H [km] = 825 T [K]/m[amu] g [cm sec -2 ]) is the atmo- 
spheric scale height [40]. Equation (2.3) represent the barometric formula 
or hydrostatic distribution in its simplest form. The barometric formula is 
also a direct consequence of a Maxwellian velocity distribution, including the 
potential energy in a planetary gravitational field. 

The total content in a column of unit cross-section of an atmosphere with 
constant scale height is given by 



poo 

= / ndz = n 0 H 

Jo 



Po 

mg 0 ' 



(2.4a) 
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The total mass of a planetary atmosphere can be expressed by 



M atm = (J j 4 t xRl (2.4b) 

where the subscript “s” denotes the values of p and g at the planetary surface 
and i?o is the planetary radius. 

A more generally applicable form of the barometric law is given by the 
differential equation 



dp d n d T dz 

7 = 7 + ^ = “7 (2 ' 5) 

which can be used throughout the barosphere. It is often desirable to take into 
account the altitude variation of the acceleration of gravity in the barometric 
law, by defining a potential or reduced altitude z’ . As long as the centrifugal 
force due to the planet’s rotation is negligible, the acceleration of gravity is 
given by 



9(z) 



9of% 

(Ro + zf 



( 2 . 6 ) 



where Ro is the planetary radius and go is the acceleration of gravity 
at the planet’s surface (go = GM/Rq , with M the planetary mass and 
G = 6.6695 x 10 8 cm 3 g 1 sec 2 ). Numerical values of go for the planets 
and relevant satellites are listed in Table 2.1. The reduced or potential alti- 
tude can then be defined by 



z' = f ^-dz = Rl f 
Jo 9o Jo 



dz 



Roz 



(Ro d- z) 2 Ro 4- z 



(2.7) 



By using the reduced altitude, the variation of g(z) is taken into account 
and go can be used in the expression for the scale height, making H truly 
constant if m,T = const. Taking into account (2.5) the atmospheric scale 
height is defined by 



H = — 



dlnp 

dz 



(2.8) 



which is the inverse logarithmic decrement of pressure. 

In the heterosphere and thermosphere the scale height varies with altitude 
as the result of the variation of m and T, leading to 



1 d H _ _1 dT _ 1_ dm _ Id g 
H dz T dz m dz g dz * 



(2.9) 



In this case a distinction must be made between atmospheric (pressure) scale 
height H and density scale height H p , since H > H p [40]. (The two scale 
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Table 2.1. Surface gravitational acceleration of solar system bodies. 



Planet 


go [cm/sec 2 ] 


Mercury 


370 


Venus 


887 


Earth 


980 


Mars 


371 


Jupiter 


2312 


Saturn 


896 


Uranus 


869 


Neptune 


1100 


Pluto 


81 


Satellites 


Io 


181 


Europa 


130 


Ganymede 


142 


Callisto 


125 


Titan 


135 


Triton 


78 


Charon 
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heights are the same only for an isothermal atmosphere with m = const.) 
This fact is important since the deceleration of artificial satellites resulting 
from atmospheric drag depends on the product pH p ' . The pressure and 
density scale heights are related by virtue of 



H P = - 



1 dp 

L p dz 



H 



1 + f3 - 2H/(Ro + z) 



( 2 . 10 ) 



where (3 = dH/dz; the term 2H / (Rq + z) arises from the variation of g with 
height. The scale heights in units of, “reduced” altitude are thus given by 



rr, H ' 

p 1 +0' 



( 2 . 11 ) 



where ft = dH'/dz'. 

Often the scale height gradient 0 can be assumed to be constant over 
certain height intervals. In this case the pressure and density distributions 
can be expressed by [41] 



( H \ 

P = Po\ir = po exp(—() 

\Ho / 

/ rr \ ^— 1 )/ (3 

n = n °[ir ) = n 0 exp[-(l + /3)C] (2.12) 

\Lio / 
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as the result of introducing a height variable £ defined by H = Ho exp(/?£) or 
equivalently d£ = dz/H, if we take (3 = dH/dz = const., so that H = Hq+/3z. 

Any of the distributions given by (2.3) and (2.12) are called hydrostatic 
distributions. They are representative of a mixing distribution when the 
scale height H corresponds to a mean molecular mass m = const. When 
each individual constituent is distributed according to its own scale height 
Hj = kBT/mjgj such a distribution is also called a diffusive equilibrium dis- 
tribution. In this case the various constituents obey Dalton’s law of partial 
pressures. 



Dynamic (Transport) Effects 



The principal transport process in the upper atmosphere is diffusion. It is due 
to a gradient in relative concentration arising from slight deviations from a 
Maxwellian distribution. The relevant properties can be described by gas- 
kinetic methods of transport theory. 

In considering the behavior of a minor constituent under diffusion, it can 
be assumed that the total pressure remains practically constant, while the 
minor constituent diffuses upward or downward due to the presence of local 
sources or sinks. This requires that the divergence of the flux Fj(z) of the 
j - th constituent obeys the continuity equation 



n r ■ 

-£=q j {z)-L j {z) (2.13) 

where qj and Lj are the sources and sinks, respectively, for the j - th con- 
stituent. 

The flux Fj which can be supported by a given density distribution rij(z) 
is given by [40, 42] 



Fj 



= n 3 W 3 = ~ n 3 D 3 



= -UjKjy 



1 d rij 1 (1 + aj) d T 

rij dz Hj T d z 

'JL dr^ JL_ 1 dT' 

nj ~dz + H + T dz . 



(2.14) 



where Wj is the flow velocity, Dj is the molecular diffusion coefficient, Hj is 
the scale height of the j-th constituent, H is the (average) atmospheric scale 
height, OLj is the thermal diffusion factor and Kb is the eddy diffusion or 
mixing coefficient. Equation (2.14) contains the effects of molecular diffusion 
(first term) as well as eddy diffusion (second term) which is of principal 
importance below the homopause, whereas molecular diffusion is predominant 
above the homopause. Figure 2.1 shows the effects of Eddy and molecular 
diffusion. 

The molecular diffusion coefficient in a multiconstituent gas is given by 
[43] 
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Fig. 2.1. Pioneer Venus measurements of He and CO 2 number densities obtained 
by a neutral mass spectrometer. The gradient of He densities above about 200 km 
altitude indicate exospheric temperatures of about 280 K. The curvature of the 
He profile below 200 km has been used after [44] to calculate the Kr> profile in 
the lower atmosphere. One can identify the homopause from altitude variations of 
both species (mco 2 = llwiHe)- According to (2.14) the homopause is defined where 
D = Kd leading to a scale height of the minor constituent at the homopause of 
H* = 2&BT/(raco 2 + rriHe)p; he. an effective mass m* = (mco 2 + uiHe)/2 [45]. 



1 _ Uk_ 

D j bjk 



(2.15) 



where the bj k are given by the binary encounters between particles of species 
j and k. Assuming rigid elastic spheres of diameters d 



_ 3 / k&T{mj + m k ) \ 2 

8 d 2 j \ 27Trrijmk ) 



(2.16) 



where d jk = \{dj + d k ). 

Thus, the molecular diffusion coefficient is proportional to 

Dj oc T 1/2 n~ l . (2.17) 

The thermal diffusion coefficient can also be estimated by assuming rigid 
elastic spheres. In the case of a minor constituent, aj -» —5/13. However, 
the rigid sphere model seems to overestimate o^; a practical value for light 
minor constituents (H, D, He) is aj — —0.25 [44]. 

There is at present no completely satisfactory physical derivation of the 
eddy diffusion coefficient K&. However, if it is assumed that K d arises from 
turbulence produced by internal atmospheric gravity waves then 
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K d = 1.4 x + X 2 z )~i 

where A x and A z are the horizontal and vertical wavelength respectively. In 
contrast to the molecular diffusion coefficient Dj ~ n -1 , the eddy diffusion 
coefficient K& ~ n~ 1 / 2 . For the case of a constant scale height gradient /?, 
the diffusive flux above the homopause can be expressed by 



Ei 

H 



dz 



+ 






Tin 



(2.18) 



Integration of (2.18) with Fj = const, leads to the density distribution under 
flow conditions [42] 
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(2.19) 



where 



rrii 



m 



+ /3otj + f3 B 



FjH 0 

n j0 Do 



C = 1 



rrij n (3 
~^ +aj(}+ 2 



The parameters in B are taken at the lower reference level (subscript 0) of the 
diffusion-controlled region, i.e., above the homopause. The bracketed terms 
in (2.19) represent a correction factor to the non-flow (diffusive equilibrium) 
condition. The diffusive equilibrium distribution corresponds to the case Wj = 
0 (cf. [41]— [52]). A more generally valid form of (2.19) for a variable flux Fj(z) 
is given by [47] 



nj (z) 




Fj(z) 

Dj ( z )rij (z) 



dz 



( 2 . 20 ) 



where rij(z) = ( To/T)rijo exp (— fj d z/Hj) is the zero flow or diffusive equi- 
librium distribution. 

A good approximation for the density under flow conditions, rij can be 
obtained in terms of the integrated column density (total content) Tt of the 
main atmosphere (density n) through which the constituent j diffuses [48]. 
The flux Fj can be expressed in terms of the ratio of density under flow to 
the zero flow (static) density r)j = rij /nj according to 






b 3 n 3 \ 

n J dz 



where bj = Djji is the diffusion coefficient at unit concentration, with n the 
density of the main atmosphere. Using n = — dVl/dz we can write 

drij Fj dVX FjVX 

a=t=^- V3 = r= + v oo 
dz bjTij dz bjfij 
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where rjoo is the ratio of actual to static density at great heights and the 
density under flow follows as 

_ FM _ 

n 3 = VjUj = -J— + VocUj . 

°3 

The flow which can be supported by diffusion in a given atmosphere without 
external sources or sinks (V • F j = 0) is limited by the condition B = C in 
(2.19) and is given by 



p* rijoDo 

Fi =_ hT 



rrij ( 1 

1 P ( OLj -F - 

m V 2 



( 2 . 21 ) 



This leads to a density distribution at the reference level (z = 0), d(ln nj)/dz = 
— (l-b/?/2 )/H. This characteristic slope is approached at lower levels for near 
maximum flow, e.g., 0.95 F? and therefore the transition from a mixing dis- 
tribution, 



d(ln rij) _ (1 + /?) 

dz ~ 

is almost indistinguishable. This is illustrated in Fig. 2.2 for hydrogen in N 2 . 
For a constant scale height it can easily be shown that the maximum upward 
flux 



n 0 A) r. nij-i 

leads to a density distribution rij = njo exp(— z/H), i.e., a distribution corre- 
sponding to the main mixed atmosphere through which the minor constituent 
diffuses. The density distribution of a minor constituent under flow and dif- 
fusive equilibrium for Hj = const, is illustrated in Fig. 2.3. The time con- 
stant for (molecular) diffusion can be derived from the continuity equation 
drij/dt = — dFj/dz , noting that r = ((l/n)(<9 n/<9£)) -1 , i.e., the characteristic 
time for attaining a diffusive equilibrium distribution is given by [42] 

T Dj ~ — = J~ n ( z ) ■ (2.22a) 

For a constant scale height, td decreases exponentially with altitude since it 
is proportional to n(z). 

Similarly, the time constant for turbulent mixing (eddy diffusion) is given 
by 



tk d ^ ^ (2.22b) 

where K& is the eddy diffusion coefficient. (The altitude where td = tk d is 
usually identified with the homopause.) 
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Fig. 2.2. Density distribution for mixing, diffusive equilibrium and flow condi- 
tions of the minor constituent hydrogen in a N 2 atmosphere, having a scale height 
gradient f3 — 0.2, including the effect of thermal diffusion (after P. Mange [42]). 




Fig. 2.3. Idealized density distribution of a minor constituent in an isothermal 
atmosphere under mixing, diffusive equilibrium and flow conditions. 
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The equation of motion of the neutral gas (such as a wind system set up 
by pressure gradients) can be expressed by [49] 

r\ i 

+ — J'inCVn - Vi) = g Vp n + 2(v n X Q) + — V 2 v n (2.23) 

Pn Pn Pn 

where v n is the (vector) velocity of the neutral gas, V[ is the velocity of the ion 
gas (which may be constrained by a planetary magnetic field), i/ in is the ion- 
neutral collision frequency, p n ,(i) is the neutral (ion) mass density, v n x fl is 
the Coriolis acceleration, with Q the angular velocity of the planet’s rotation 
and 77 the coefficient of viscosity (r]/p n is the kinematic viscosity). The second 
term on the left-hand side is referred to as the ion drag term, which arises from 
the interaction between the neutral atmosphere and the ionosphere. Attempts 
have been made to deduce the wind-system in the terrestrial thermosphere 
arising from the diurnal pressure field [49] . 

Figure 2.4 shows the thermospheric wind field derived from (2.23). The at- 
mospheric motion (winds) must also satisfy the conservation of mass through 
the equation of continuity 



dp 

dt 



V • (pv) 



for the neutrals and simultaneously for the ions. Accordingly, vertical motions 
must occur which balance the divergence or convergence resulting from the 
horizontal wind field. 

Atmospheric disturbances can propagate from the lower to the upper 
atmosphere where they may be responsible for oscillations in the neutral 
gas as well as in the ionospheric plasma in form of travelling ionospheric 
disturbances (TID) [50]. 

Oscillatory solutions can be obtained from the equations of motion and 
continuity for the neutral gas in (2.23), having the form 



Ap Ap 

p p 



ocv n (x exp i(t at — kr) 



(2.24a) 



where k = (k x + k z ) 2 is the complex wave number, with k x , z = 27t/A x?2 . For 
upward propagation in a lossless atmosphere one obtains the solution 



Ap 



(x exp 




(2.24b) 



since the atmospheric scale height H and the sound speed c s = ( 7 k&T/m)% 
are related by H = cl/'yg, where 7 is the ratio of the specific heats. This 
result is a direct consequence of the condition that the energy flux remain 
constant, i.e., 



\pv\ = const. 
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Global wind field 
Height 300km Equinox 




Fig. 2.4. Simplified global thermospheric wind system at equinox derived from 
(2.23) using model values of Vp n /p n from satellite drag data and ion drag based on 
a semi-empirical model ionosphere (maximum wind speed ~ 300 m/sec) (courtesy 
of P. Blum and I. Harris). 



since in an isothermal atmosphere p oc exp(— z/H). Waves for which such 
propagation can take place are determined by the following cutoff condi- 
tions. Waves whose angular frequency lj > u g = 7 g/ 2 c s are termed acoustic 
waves , whereas those for which uo < 00 g = (7 — 1)2 g/c s are termed inter- 
nal atmospheric gravity waves. The cutoff frequency o; a is called the acoustic 
low frequency cutoff ; while the high frequency cutoff for internal atmospheric 
gravity waves u g is called the (isothermal) Brunt-Vaisala frequency. Internal 
atmospheric gravity waves can be considered a type of low frequency acoustic 
wave. (Tidal waves are essentially internal atmospheric gravity waves where 
the effect of the Coriolis force is included and which propagate primarily hori- 
zontally.) These waves can be produced from energy supplied by tidal forces, 
from large-scale (planetary) wind systems, or as the result of atmospheric 
heating by energetic particles as occurs in the auroral zone. 

In reality, upward propagation of internal atmospheric gravity waves for 
a given uj is limited by the kinematic viscosity of the atmosphere as well as 
by reflections from a thermal barrier (i.e., where dT/dz > 0) so that only 
certain wavelengths are able to penetrate to higher altitudes [51]. 
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II.2 Thermosphere: Thermal Structure 

The energy budget of the upper atmosphere of planets is primarily governed 
by the heating of the gas due to the absorption of solar extreme ultraviolet and 
X-ray (XUV) radiation by atmospheric constituents, by heat transport due to 
conduction and convection and by heat loss due to emissions in the infrared 
(IR) [52, 53]. Radiative loss by IR occurs when atmospheric constituents are 
present which have transition levels in the infrared. This is the case for atomic 
oxygen in the terrestrial atmosphere and H 2 , CO and CO 2 in other planetary 
atmospheres. 

In addition to the XUV flux there are other possible sources of heating of 
the upper atmosphere of planets. Among them are collisions between charged 
particles and neutral constituents, Joule heating, conversion of dynamic en- 
ergy into heat, such as from tidal motions and internal atmospheric grav- 
ity waves, absorption of hydromagnet ic waves and either direct or indirect 
heating by the solar wind. Direct solar wind heating can occur in planetary 
atmospheres without a significant screening magnetic field (magnetosphere) 
as on Venus or Mars, while indirect effects can occur via magnetospheric 
processes which are again induced by the solar wind (e.g., heating associated 
with auroras). 

The thermal balance in the thermosphere is governed by a heat conduction 
equation of the form [49] 

r dT 1 

pc v -+v n -VT + pV-t; n - V • (K n VT) = Q xvv - L m (2.25) 

where p is the atmospheric mass density, c v is the corresponding specific heat 
at constant volume, v n is the velocity and p the pressure of the neutral at- 
mosphere, K n is the thermal conductivity which is a function of temperature 
according to K n = AT 0 T S , Qxuv represents the total volume heat production 
rate, which may be the sum of a number of processes, including recombi- 
nation heating and heating by collisions with charged particles of the form 
Q [n = A(T e — T n ) + B(T[ — T n ) where A and B depend on the appropriate 
collision frequencies between electrons and neutrals and ions and neutrals, 
respectively. Lir, the radiative loss term, depends on the number density of 
the emitting state which is related to the total number of the emitting species 
n x by a partition function f(T) and is thus of the form [52] 

Lir = n x /(T); 

the details have to be derived from radiative transfer considerations [53]. 

A simplification of (2.25) has been considered by a number of workers in 
some detail, and has the form [54] 

(dT dT\ d / frri .dT\ „ 

^ af + *"■ & ) - Tz r " (T) aJ ) = Qxw “ Lm 
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together with a diffusive equilibrium distribution for n. For the one-dimensio- 
nal case (vertical variability z only) the condition dp/dt = 0 allows replace- 
ment of c v by c p and the expression of v nz as the velocity of a parcel of air 
identical to the vertical velocity of the isobaric surface ( “breathing velocity” 
of the atmosphere) given by 



^ N f z 1 dT , 

Actually, v nz must be considered the sum of the breathing velocity and a 
divergence velocity arising from the horizontal wind system by virtue of con- 
tinuity [55]. 

The principal heat source for the thermosphere is generally considered 
to be the conversion of absorbed solar XUV radiation into thermal energy 
by superelastic collisions. The heat production rate due to the absorption of 
solar XUV radiation for a given wavelength and constituent is given by 

Qxuv = e j n j aJxvv^~ T (2.26) 

where £j is the fraction of absorbed XUV energy which is transformed into 
thermal energy, i.e., the photoionization heating efficiency , which is thought 
to lie in the range from 0.3 to 0.6 and may also differ for different spectral 
ranges; lx uv is the energy flux [erg cm -2 sec -1 ] at the XUV wavelength range 
outside the atmosphere and is related to the photon flux $oo by 




where f) is Planck’s constant, c the velocity of light; <r a is the appropriate 
absorption cross-section (typically in the range 10 -18 -10 -16 cm 2 ) and r is 
the optical depth defined by 

poo 

r = / sec% • 0 a • Uj dz (2.27) 

with % the solar zenith angle (cosx = sin <p sin 5 + cos(pcos<Scos(i?£), where 
ip is the latitude, 6 is the solar declination, ft is the atmospheric (planetary) 
rotation rate and t is time measured from the noon meridian); for x > 75°, 
secx must be replaced by the Chapman function Ch(x) which allows for the 
curvature of the planet (cf. Sect. LI). The solution of the time-dependent 
heat balance equation (2.25) leads to a time variation of the temperature in 
the thermosphere, having a maximum which occurs later than that of the 
XUV heat input, due to the thermal inertia of the atmosphere. Figure 2.5 
illustrates the diurnal variation of the exospheric temperature and that 
of the terms in the heat balance equation represented by their integrated 
values, i.e., the total XUV heat input, the total IR loss, the total conduction 
loss and the mechanical energy f pV • v n dz, as obtained from a numerical 
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THERMOSPHERIC ENERGETICS: 




Fig. 2.5. Diurnal variation of Earth’s exospheric temperature Too and the inte- 
grated values (above 120 km) of the terms in the heat balance equation. “Ionization” 
due to XUV radiation contributes to thermospheric heating. “Photodissociation 
energy” represents the integrated heat input due to dissociation by the Schumann- 
Runge continuum; “Photoelectron energy” represents the heat input due to XUV 
photoelectrons; “Plasma feed-back” includes the heat source due to collisions be- 
tween the neutral and electron-ion gas and due to exothermic chemical reactions 
(dissociative recombination); “Mechanical energy” represents the integrated heat 
input due to the expansion and contraction (f pV • v n d z) of the atmosphere; “Ra- 
diation loss” represents the IR radiation loss above 120 km and “Conduction loss” 
represents the integral of the heat conduction term (courtesy of S. Chandra). 
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solution of (2.25). In addition to the XUV heat input, collisions between the 
neutral and the electron-ion gas can provide an additional heat source for 
the thermosphere according to 



Qcoii = 2iVi 



A*in l 
‘ mim n 



^k B {Ti - T n ) + ^m n (v ; - v n 



as can exothermic chemical reactions such as dissociative recombination (cf. 
Chap. IV) 



XY+ 4- e V + Y + AE. 

Figure 2.6 shows the vertical temperature profile at the diurnal maximum 
and minimum for the conditions shown in Fig. 2.5. A problem has arisen 
for the terrestrial thermosphere which has been termed the diurnal phase 
anomaly between atmospheric density and temperature. From satellite drag 
observation it was established [56] that the atmospheric density in the ther- 
mosphere reaches a maximum around 1400 LT (local time) and consequently 
a temperature maximum was inferred for the same time from atmospheric 
model considerations. 

The magnitude of the diurnal phase anomaly is a function of atmospheric 
composition, since the individual atmospheric constituents attain their di- 
urnal maxima at different times, e.g., He in the morning hours and O in 
the early afternoon. Thus, the lighter the constituent, the larger the phase 
anomaly. The solution of the heat conduction equation (2.25) on the other 
hand leads to a maximum in the temperature at about 1700 LT. A number 
of suggestions have been made to explain this phase discrepancy which range 
from the original assumption of a secondary ad hoc heat source [54] to the 




TEMPERATURE (K) 



Fig. 2.6. Vertical temperature profiles in the terrestrial thermosphere at the time 
of diurnal minimum and maximum. 
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Table 2.2. Dependence of diurnal variation of exospheric temperature Too on basic 
input parameters. 



Cause 
Increase of: 


Effect 

T max 


T min 

To 


fmax 

T min 


Rotation rate 

of atmosphere, Q r 


(superrotation) 


-* 


T 


rv 


Specific heat c p 


Thermal conductivity K n 


e- 




rv 


Absorption cross section, a a 


<- 


t 


t 


Pressure at base level, po 




t 


rv 


Solar activity, loo 


$ 


t 


$ 



t increase, ^ decrease, rx goes through max, $ ~ const, earlier, — * later. 



inclusion of dynamical effects and extension of the one-dimensional heat con- 
duction to two and three dimensions [57, 58]. In addition, the ratio between 
the maximum and minimum temperature also differs between theory and that 
derived from satellite drag observations, the latter being 1.3, whereas the for- 
mer is ~ 1.5-2. Thermospheric temperature data inferred from ground based 
incoherent scatter radar measurements suggest the peak temperature occurs 
near 1700LT with a diurnal amplitude closer to the theoretical value. In the 
light of this result, a suggestion has been made [59] that the phase anomaly 
is due to the fact that the density and temperature at the homopause is not 
time independent but show a diurnal variation, whereas others [60] view the 
phase anomaly as the result of wind-induced diffusion. A parametric study 
of the thermospheric heat balance equation of the type 



pc p 



(fT_ 



dT 



dT 



dt +VB dz 



— Q + 



dz 



dT 

n dz 



(2.28) 



suggests that the diurnal variation of responds to changes in the basic 
input parameters of (2.28) as shown in Table 2.2 [61]. One should also note 
that super-rotation of planetary atmospheres may occur; i.e., the atmosphere 
(or parts of it) may rotate faster than the planet [62]. 

When a two-dimensional atmospheric model including winds is considered 
[61], the “diurnal amplitude” of temperature is found to be 



Tm ax Q(f2) 

^min |-^con — fiC\ 

where Q(Q) is the diurnal XUV heat source and Q = 27r/d the angular fre- 
quency, with D the solar day, i.e., in general, i? is the rotation rate of the 
planetary atmosphere; L con is the heat loss due to conduction and convection, 
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i.e., mass transport by winds, both assumed to be proportional to temper- 
ature and C = c v p is the heat capacity of the atmosphere. The time of the 
diurnal temperature maximum is related to these quantities by 

t 

‘'max ^ j 

-t>con 

Empirical relations have been derived for the variation of exospheric tem- 
perature with solar activity and geomagnetic activity. As an indicator of the 
solar XUV intensity, the 10.7 cm solar radio flux (2800 MHz) observed on 
Earth has generally been employed, since continuous observations over long 
time intervals of the solar XUV radiation are very difficult. However, the 
close correlation between solar XUV and the 10.7 cm radio flux, F 10 . 7 , has 
been established over limited time intervals [63]. Typical values of F 10.7 in 
units of 10 -22 Wm -2 Hz ” 1 averaged over several solar rotations are: 75 for 
low solar activity (solar minimum), 150 for medium activity and 230 for high 
solar activity (solar maximum). 

Empirical relations between the terrestrial exospheric temperature de- 
rived from satellite drag and the 10.7 cm flux have been developed from 
observations extending over a full 11-year solar cycle [56]; a nearly linear 
correlation was found between T ^ and F 10.7 (averaged over several solar 
rotations = 27 days), according to which 

AT oo _ 4K 

AF10.7 ~ 10 -22 W m ~ 2 Hz -1 

The temperature change due to magnetic activity is 

ATqq = 20 K p + 0.03 exp K p 

where K p is the planetary magnetic index, which in turn is related to the 
solar wind velocity u sw according to 

v sw [km/sec] = 330 + 67.5 K p . 

It should be noted that for planets without an appreciable intrinsic magnetic 
field, more direct solar wind heating may occur. 

More recently it has been recognized that E 10.7 is not a perfect proxy 
for the XUV radiation responsible for heating and ionization of the upper 
atmosphere. Particularly at high solar activity E 10.7 seems to overestimate 
the XUV flux. Studies of the solar control of the ionospheres of Venus [64] and 
Mars[65] suggest that lx uv ^ 10.7 with m ~ 0.9; i.e. a slightly non-linear 

dependence. A comparison between T 10.7 as proxy for the solar XUV input 
on the terrestrial thermosphere has led to a corrected E 10.7 flux, called E 10 . 7 , 
as proxy for the integrated energy flux in the XUV range [71]. 

In the following we shall illustrate some of the general characteristics of 
planetary thermospheres which are the result of solar XUV heat input and 
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heat conduction. For this purpose we consider a simple form of the heat 
balance equation neglecting the IR heat loss. For steady state conditions 
(■ dT/dt = 0) we have 



d z 



^(* n S> =Qxuv 



and by integration 

/»oo 

S = / Qxuv dz = —CX£jIx uv(l - e ~ T ) 



noting from the definition of optical depth that dz = dr /nja a . Since the heat 
flux is S' = —K n j(T) dT/dz, we can write 



d T _ Sjl x uv(l — e r ) 
dz ~ aK n j(T) 



(2.29) 



This represents the thermospheric temperature gradient which can be sup- 
ported by the XUV heat input and heat conduction. 

The factor a takes on values appropriate for the rotation of the planet: 
about 4 for rapidly rotating planets like Jupiter and about 2 for slowly rotat- 
ing or tidally locked planets, because of the ratio to the cross-sectional area 
of the surface of the planetary sphere. The value of a may be reduced close 
to 1/4 by the action of strong thermospheric winds expected on exoplanets 
at close orbital distances. 

In the exosphere, where the density is low, r — > 0 which leads to dT/dz = 
0, i.e., an isothermal atmosphere. This constant temperature is also called 
the exospheric temperature 



T{z) = Too{l - aexp(— 6z)} 



where 



a = 



Too — To 



and 



b = 



1 



~ Tq 



d T 



dz 



Zo 



The temperature gradient will have its maximum value where the optical 
depth is large (r — oo), i.e., in the lower thermosphere. This maximum 
temperature gradient is given by 



/ dT\ £jl x uv 

V^/ma x ~Knj(T)- 
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Table 2.3. Solar energy flux of a 4.5 Gyr old G-type star (Sun) outside a planetary 
atmosphere normalized to that observed in the terrestrial exosphere as a function 
of orbital distance. 



Planet 


/xuv(P) 

7xuv(e) 


d [AU] 


OGLE-TR-56 b 


1890.35 


0.023 


HD209458 b 


493.82 


0.045 


Mercury 


6.65 


0.387 


Venus 


1.9 


0.723 


Earth 


1.0 


1.0 


Mars 


0.43 


1.52 


Jupiter 


3.7 x 10“ 2 


5.2 


Saturn 


1.1 x 10“ 2 


9.54 


Uranus 


2.7 x 10 -3 


19.18 


Neptune 


1.1 x 10 -3 


30.06 


Pluto 


0.65 x 10“ 3 


29.4 



The thermospheric temperature gradient is therefore directly dependent on 
the XUV heat input (and the heat conductivity). For medium solar activity, 
the total flux at Earth in the wavelength range 1310 - 270 A (excluding the 
strong Lyman-a line, 1216 A) is about 2 erg cm -2 sec -1 (for further details 
see Chap. I). 

The energy flux for a given wavelength A outside the planet’s (P) at- 
mosphere can be scaled from the terrestrial values according to I^(P) = 
7^(®)/d 2 , where d is the mean distance of the planet from the Sun in astro- 
nomical units (AU). Table 2.3 shows XUV energy fluxes at various objects 
and orbital distances. The time constant for equalization of a temperature 
difference by means of heat conduction is given by 



Tcond j 



njC v k^L 2 

Knj(T) 



(2.30) 



where L is a scale length over which this temperature difference exists. 

By performing a second integration of the simple heat balance equation 
we can also determine the variation of the “exospheric” temperature with 
the variation in solar XUV heat input over a solar cycle [67]. 

The integrated form of (2.29) is given by 




— e r )d z 



where is the exospheric temperature, To is the temperature at the 
mesopause (which depends strongly on the IR loss) and z and zq are the 
corresponding altitudes. If we assume that r -4 0 at i.e., in the exo- 
sphere, and t — y oo at zq, i.e., at the base of the thermosphere, and allowing 
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Table 2.4. K n j = KqjT s ^ [ergs cm - 2 sec -1 °K _1 ] 



Constituent 


K 0j 


Sj 


n 2 


27.21 


0.80 


o 2 


18.64 


0.84 


C0 2 


1.5 


1.23 


CO 


26.16 


0.80 


0 


67.1 


0.71 


He 


20.92 


0.75 


H 


16.36 


0.73 


Ar 


18.3 


0.80 



for (nt^oojaa) 1 z — zq (since the absorption cross-section cr a ~ 10 17 cm 2 
and exospheric densities are of the order of 10 8 cm -3 ) we obtain 



Ko (IS 



-Tn 



Tlj (^oo)^a 



For rij{z 00 ) we can substitute the density at the base of the exosphere n c = 
(crif) -1 (see next section), where H = k&Too/mg is the exospheric scale 
height, and a is the gas-kinetic collision cross-section. 

The variation of Too with solar XUV energy flux is therefore determined 
by the temperature dependence of the thermal conductivity. It has been 
generally assumed, based on a rigid sphere approximation, that K n j = AjT 5 . 
(Note that the heat conductivity is related to viscosity 77 by K n = |7/c v , cf. 
[43].) However, experimental data [68] and quantum mechanical calculations 
[69] show that for atmospheric gases 



K nj = KojT Sj 

where Sj > The heat conductivity and its temperature dependence for 
likely constituents of planetary atmospheres is given in Table 2.4. 

From the effective heat production below the exosphere, which is balanced 
by the divergence of the conductive heat flux of the XUV radiation, one gets 
the following expression for [70] 



rpsj _ Sj £j fee 

00 K 0j a rrij g 



{ E [r(A)] + In [t(A)] + 7 — fl 
l m k L 



r* 2 

/ ^xuv(A) x 
J Ai 

1 dA + lS', 



-r(A) 



(2.31) 



where A is the wavelength, r is the optical thickness, E is the exponential 
integral, 7 is Euler’s constant, rrij and m k are the masses of the atmospheric 
constituents and To is the temperature at the base of the thermosphere. In 
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the case of a hydrogen dominated thermosphere j = k. After some approxi- 
mations, one obtains 



£j lx uv k B a Cj , ^ 

Ty ~ 0 

a K 0j rrij g a &j 



(2.32) 



in which a Cj and cr aj are the collision and absorption cross-section, respec- 
tively. Tq is negligible in cases where T^ is large. If planets 1 and 2 have an 
identical gas composition j in their thermosphere one obtains the following 
scaling relation [71, 72] 



(Toi - Tp ) 1 _ /xUViff2 
(T%o ~Tq) 2 Ixu v 2 9i 



One can see from 2.33 that T^ depends on lx uv? which decreases with 
distance from the star, and on g , which is related to the mass and size of 
the planet. Thus, 2.33 allows us to estimate T^ and the thermal escape of 
hydrogen-rich exoplanets as a function of orbital distance, planetary mass 
and planetary size. 

From satellite drag observations it is known that the terrestrial exospheric 
(daytime) temperature ranges from ~700 K at solar minimum to ~2000 K at 
solar maximum, i.e., a variation by a factor of about 2.5 while the exact varia- 
tion of the XUV energy flux over a solar cycle has yet to be established. Since 
atomic oxygen is the principal constituent of the terrestrial thermosphere, ac- 
cording to oc lU°' 71 one can infer from the temperature ratio of 2.5, that 
the XUV energy flux should have changed by a factor of ~ 1.9. From studies of 
the solar control of the ionospheres of Venus and Mars [64, 71] the dependence 
of neutral scale height H at ionospheric levels on solar activity expressed by 
Tio .7 for the solar XUV input was found to be H oc F [ g 7 with m « 0.15 for 
both planets. Since the scale height H = k&T/mg , this dependence can be 
thought to hold for T^ oc T{q 7 . On the other hand according to the sim- 
plification of (2.32) T^ oc (Ix\jv/rng) 1/s - The mean molecular mass m in a 
multicomponent thermosphere is subject to diffusion and will vary with tem- 
perature and thus with solar activity, increasing at the thermospheric levels 
according to m oc n b.7- Since lx uv oc F®q 9 7 one obtains for the variability of 
exospheric temperature with the F 10.7 proxy 



-,( 0 . 9 -k)/s 



where from the exponent m « 0.15 one finds for Venus and Mars (s = 1.23) 
the variability of the exponential of the mean molecular mass according to 
k & 0.7. For Earth the dependence of T^ on T 10.7 was empirically found from 
the analysis of a series of satellite missions to be T^ oc F ^ 2 7 . By making the 
same comparison of exponents as above one finds that for s = 0.71 a value 
of k « 0.7 is found, i.e., a similar dependence as on Venus and Mars. 
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The exosphere temperature T ^ is the characteristic planetary tempera- 
ture for thermal escape and is thus an important factor for the evolution of 
planetary atmospheres and their water inventories. In contrast to the surface 
temperature, it is not directly correlated to the distance from the radiation 
source because of the dependence on both energy input and heat conduction. 

The effective temperature T e ff is generally the result of equilibrium be- 
tween short wave radiation according to the source absorbed by the surface 
and emission of the long wave radiation according to the Stefan-Boltzmann 
law. For an energy flux 5 at a distance d from the source (Sun, star, etc.) 



5( 1 - A)nR 2 0 
d 2 






and thus 



T e ff = 



' 5(1- A) ]* 
aasnd 2 



where A is the planetary albedo and ctsb the Stefan-Boltzmann constant. 
This radiation surface temperature is only equivalent to the actual surface 
temperature T s when the atmosphere of a planet is transparent to both the 
incoming short wave radiation (visible) as well as the outgoing infrared ra- 
diation (IR) in consequence of Wien’s displacement law. If the atmosphere 
contains polyatomic gases capable of absorption and emission in the IR, the 
greenhouse effect will ensure making the actual surface temperature higher 
than T e ff 



Ts = Teff (1 + Tir ) 4 



where tir is the effective IR optical depth due to the presence of greenhouse 
gases. 

The consequences of high exospheric temperatures for the escape (expan- 
sion) of gases from planetary atmospheres will be discussed in the following 
section. The escape becomes excessively large for a critical temperature 



T C rit — 



3gm /0 

w ( 0+ } 



corresponding to the condition when the thermal energy of the gas kinetic 
motion (3&bT)/2 overcomes the gravitational energy corresponding to an 
escape parameter X = ( Voo/vo ) 2 = 1.5 where Vqq is the escape velocity 
and vo is the most probable velocity of a Maxwellian velocity distribution of 
temperature T cr i t . This condition was also called blow off by the astrophysicist 
E. Opik, who first recognized the importance of hydrodynamic escape [73]. 
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II. 3 Exosphere: Atmospheric Escape (Evaporation) 

At any level of a planetary atmosphere there will be some particles which 
move upward with a velocity greater than that required for escape from the 
gravitational attraction of the planet. These particles are representative of the 
high velocity tail of the Maxwellian velocity distribution. The actual escape, 
however, will depend on how frequent (or rather how infrequent) collisions 
are with other particles. At high enough altitudes the collision frequency will 
be low enough (i.e., the mean free path long enough) for escape to take place. 
This region is called the exosphere. The probability for a molecule travelling 
a distance z without a collision is given by 

P(z) = exp (-|) , 

where A is the mean free path 

A = (ncr)” 1 

with n the total number density and a the gas-kinetic collision cross-section. 
Thus, the mean free path represents the distance for which this probability is 
1/e. The base of the exosphere z c (exobase), originally called the critical level, 
is usually (but somewhat arbitrarily) defined by this probability requiring 
that 

n(z)a dz = n c Ha = — = 1 . 

A 

Exospheric conditions, i.e., the possibility of evaporation, are therefore as- 
sumed to hold at altitudes where the mean free path is equal to or greater 
than the local scale height, A > H. The density at the exobase is given by 
n c = (aH)- 1 . This density value then defines the critical level or the exobase. 
Numerically, n c [cm -3 ] « 3 x 10 9 /H [km] for a a « 3 x 10 -15 cm -2 and the 
total content of an exosphere is 

Woo = <J ~ 1 ~ 3 x 10 14 cm -2 . 

A “sharp” exobase obviously represents a highly idealized concept. Accord- 
ingly, it is assumed that below and at the critical level a Maxwellian velocity 
distribution prevails, whereas in the exosphere the velocity distribution is 
truncated, due to absence of particles with speed above the escape velocity. 
Thus, the barometric law applies below the exobase, i.e., in the barosphere, 
but does not hold strictly in the exosphere. For this reason, the exobase may 
also be called the baropause (the pressure at this level is « mg^l <*;). 

The minimum velocity for escape Voo is that for which the kinetic energy 
of the particle balances the potential energy in the gravitational field 

(2 GMY 

v oo — f ) — V^gR 



dz 

L m 
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Table 2.5. Escape velocities from solar system bodies. 



Planet 


Voo [km/sec] 


Mercury 


4.3 


Venus 


10.4 


Earth 


11.2 


Mars 


5.1 


Jupiter 


59.5 


Saturn 


35.5 


Uranus 


~ 21.3 


Neptune 


- 23.3 


Pluto 


1.2 


Satellites 


Io 


2.56 


Europa 


2.02 


Ganymede 


2.74 


Callisto 


2.45 


Titan 


2.64 


Triton 


3.1 


Charon 


0.75 



where G is the universal gravitational constant ( G = 6.6695 x 10 -8 cm 3 g -1 
sec -2 ), M is the mass of the planet and R is the radial distance from the 
center of the planet. 

Table 2.5 lists the escape velocities for the planets and relevant satellites 
for R = Ro; note that the escape velocity is y/2 times the circular orbital 
velocity. In the Russian literature the orbital velocity is referred to as the 
“first cosmic velocity” , while the escape velocity is called the “second cosmic 
velocity”. Our discussion of escape from planetary atmospheres is based on 
the classical Jeans treatment [74], with extensions based on other works [75, 
76, 77], primarily in conjunction with the density distribution in a planetary 
exosphere. The basic assumptions of this treatment are: a) an isothermal 
atmosphere, b) a Maxwellian velocity distribution even at the escape level. 

The outward flux of particles whose velocity exceeds the local velocity of 
escape at a planetocentric distance R, i.e., the escape flux F 0 0 (also called 
the Jeans escape flux) is obtained by the product of the Maxwell-Boltzmann 
velocity distribution function and the vertical component of velocity (v > Voo) 
over the upward hemisphere. 

The escape flux relating to the critical level i? c , (exobase) is given by the 
Jeans formula 



F — 

-*•00 — 



_vo_ 

2 0 ? 



• n ( 



CJ 



(1 + X c )e 



-x c 



(2.34) 
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where n C j is the density of the escaping constituent at the exobase, Vo = 
{2k^T 00 /nrij) 2 is the most probable velocity at a Maxwellian velocity distri- 
bution with a temperature (^0 = 1.3 x 10 4 (Too/m) 1 / 2 [cm/sec], with m 
in units of amu) and X is the escape parameter defined by 



X(R) = 



GMrrij 

Rk^Too 



(2.35) 



The exobase is defined by the total density at that level, n c = • n C j. Thus, 

X c is the escape parameter at the exobase. The escape flux at any level above 
the exobase can be expressed by 

( Xl. \ 2 

J * 



This flux decreases as R ~ 2 to satisfy continuity; (vr) is the so-called effusion 
velocity which is defined by the ratio of escape flux to the concentration of 
the escaping constituent. It should be noted that (vr) is always subsonic. 

The escape flux Foo is related to the maximum diffusive flux F* of the es- 
caping constituent in the barosphere by virtue of RfF* = RfnjWj = R^F^ = 
n C jRl(v r) where Ri is the lower boundary (Ri < R c ); i.e.; n ]C depends on the 
diffusive flux in the barosphere. If Wj < (vr), then an “escape bottleneck” 
occurs leading to an adjustment of density. 

The above expressions for the escape flux are based on the assumption of 
a Maxwellian velocity distribution at the exobase while escape itself leads to 
a perturbation of the Maxwellian velocity distribution. A statistical Monte 
Carlo treatment however shows that the simple Jeans formula in (2.34) over- 
estimates the escape flux at worst by about 30% [77, 78, 79]. As long as the 
departure from a Maxwellian velocity distribution is not too great ( X >> 1.5), 
X can also be expressed as a height parameter X = R/H. It should be noted 
that escape is important for X < 15, an atmosphere is completely bound 
to a planet for X « 30, while for X < 1.5 the exosphere becomes unstable 
and escape may become almost arbitrarily high [73]. The latter corresponds 
to the case when the thermal energy of the gas kinetic motion, (3/2)A;bT, 
exceeds the gravitational energy, i.e., u 2 /v ^ < 1.5. 

Remember that uq : u : u rms = y/2 : \J S/n : \/3; the mean square veloc- 
ity is u 2 = |uq; and u rms = (u 2 ) 2 . Typical values of the escape parameter 
X c and for the most probable velocity for the lightest constituents of plane- 
tary atmospheres are listed in Table 2.6. A useful formula for the numerical 
evaluation of the escape parameter is 



X = 60 m [amu] • v^R) • 

On a rotating planet, particles having velocities in the direction of the rota- 
tional motion will more easily attain escape velocity than those moving in 
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Table 2.6. Escape parameter X c for hydrogen. 
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10 




6 


20 


12 


Jupiter 




~1400 















Sm: Solar minimum conditions; SM: Solar maximum conditions. 



the opposite direction. The resulting escape flux will therefore be primarily in 
the forward direction; the magnitude of the escape flux at the equator is also 
several times greater than that at high latitudes. The parameter governing 
rotational effects is given by 



/' m 3 n 2 R 2 y = OR 

\ 2/CbToo ) Uo 



(2.36) 



where ft is the angular rotation velocity. 

The flux ratio between the pole and the equator for large Y is enhanced 
for increasing X. The difference between the fluxes at the pole is significant 
for different X values and is greater than that at the equator [80, 81]; the 
escape flux is always greater for smaller X (more efficient escape) . The effect 
of rotation is illustrated in Fig. 2.7. 

The concept of escape time is related to the available (total) column 
density for escape, The ratio Oloo /n C j, can be called the “escape length”, 
corresponding to the length of a square centimeter column of density n C j 
containing Ttoo particles. The loss rate of OToo is given by 

= ~ F °° = -taKj = -foOWoo • ( 2 - 37 ) 



Integration of (2.37) gives 

Woo (t) =%x>(t = 0)e-‘/ Tj 



where the Jeans escape time is 



tj 



9too 

(n cj (v r )) 



(2.37) 



Thus, Tj represents a time constant and the fraction of the total content 
lost by escape is given by (1 — e -t / Tj ). If tj to (for to == 4.5 x 10 9 yr = 
1.4 x 10 17 sec), e.g., to = with k > 7, then tj can also be considered as 
the time for (virtually) total escape, since e _t °/ rj < 10“ 3 . 
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Fig. 2.7. Effect of planetary rotation on the escape flux as a function of latitude, 
for parametric values of the escape parameter X = ( Voo/uq ) 2 and the rotational pa- 
rameter Y = QR/uq. Normalized escape flux corresponds to nj C (vR) /uo (courtesy 
of R. E. Hartle). 



Normally the big uncertainty lies in the estimate of 9too, since the total 
content available for escape may include the constituent in non-escaping form, 
e.g., for hydrogen, water in an evaporating ocean. 

Under the assumption that 9too ~ n C j • H , the escape time is given by [75] 



~ 1 / 27^00 e* c 
g \ rrij ) 1 + X c ' 



(2.38) 



Density Distribution in an Exosphere 

In computing the density distribution in an exosphere, not only the escape 
component but also the particles in elliptic ballistic orbits, i.e., those in- 
tersecting the exobase and in bound (satellite) orbits (not intersecting the 
exobase) must be considered. Because of the escape component, the velocity 
distribution in the exosphere departs from a Maxwellian velocity distribution 
and the barometric law ceases to describe the density distribution perfectly, 





II.3 Exosphere: Atmospheric Escape (Evaporation) 57 




Fig. 2.8. Group of exospheric particles in different orbits: 1) ballistic orbits v < Voo] 
2) bound satellite orbits; 3a) ballistic escape orbits and 3b) incoming capture orbits 
(v > Uoo); and 4) hyperbolic orbits (interplanetary particles) (after F. S. Johnson 
[83]). 



although it can provide in some cases a good approximation to altitudes well 
above the exobase. The contribution of the various exospheric components 
is derived by determining the fraction of a Maxwellian velocity distribution 
due to these components. The particle population in the exosphere can be 
divided into the following components as shown in Fig. 2.8: 1) ballistic orbits 
(v < Uoo); 2) satellite (bound) orbits ; 3) ballistic (escape) orbits (3a) and in- 
coming capture obits (3b) {y > Voq) and 4) hyperbolic orbits not intersecting 
the exobase (interplanetary particles) . The proportion of particles falling into 
each group can be obtained by integrating the velocity distribution function 
over appropriate ranges, represented by a partition function £ [41, 82, 83]. 
The density at any level R in the exosphere can be expressed by the product 
of the “barometric” density nMj (based on a complete Maxwellian velocity 
distribution), and the partition function £, according to 

nj{R) = n Uj {R)Z{X) = n cj e~^- x) ■ £(X) . (2.39) 

For no dynamical restriction of orbits £ = 1, nMj = / 0 °° f(v)dv. 

The complete distribution of particles in the exosphere, normalized to the 
barometric density, is given by 



£ — Ci + £2 + £3 + £4 = 1 



where the subscripts refer to the components of the total particle population 
identified above and illustrated in Fig. 2.8. Particles related to the critical 
level are represented by the partition functions [41, 76] 



£1 + £3 = 1 — 



(Xc 2 -* 2 ) 



exp 



X 2 

(X + X c ) 



} 



(2.40) 



and particles not reaching the critical level (exobase) by 
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. * {X 2 -X 2 )i 

&+U = t; exp 






X “ r 1 (X + X c 
The elliptic components ( v = v, x ) can be expressed by 

a + 6 = ${x$) + ±<P"(xi) 



(2.41) 



(2.42) 



where <£(X) is the error function and its second derivative, which are 
tabulated. The hyperbolic components (^>^ 00 ) are expressed by 



e 3 + ^ = l-#(xi)-^ , (xi). 

The bound (satellite) component is given by 

X 2 



t (X c 2 -X 2 )i 
6 = - £— 77 — exp 



X c 



{-aSnl 






X 



X 



(X + X c )i; + 2^ V(X + X C )5 



while the hyperbolic interplanetary component is given by 

X 2 



, (X 2 - X 2 )i 

U = 77 ex P 



X c 

1 -<£ 



(-B&1 



(x + x c ) 

X \ 1 



' * 

2 V (X + X c )s 



k (X + X c )* 

The above equations lead to the components at the exobase 

£1 = *(X$) + ±&'(X*) i . . . (t; < Voo ) 
6 = 1 - $(X?) - . . . (v > Voo ) 

f 2 =X 4 = 0. 



(2.43) 



(2.44) 



(2.45) 



(2.46) 



Equations (2.46) show that £3 is practically negligible for a constituent whose 
escape rate is very small, i.e., which has a large value of X c . Figure 2.9 
shows exospheric distributions of various components, represented by the 
partition functions £ given above, for hydrogen in the Earth’s exosphere at 
Too = 1500 K (X c = 4.65). The relative concentrations of terrestrial hydro- 
gen and helium at = 1700 °K are shown in Fig. 2.10 [41]. The hydrogen 
distribution based on the ballistic and escaping components departs strongly 
from the hydrostatic (barometric) law (X c = 4.35), whereas the departure 
is almost non-existent for helium at the same temperature (X c = 17.4). For 
a constituent subject to escape, the density decreases with increasing tem- 
perature, in contrast to constituents controlled by the barometric law. The 
fraction of particles with velocity v > as a function of the escape param- 
eter X is illustrated in Fig. 2.11, and Fig. 2.12 shows the escape parameter 
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Fig. 2.9. Partition functions £ for the group of particles shown in Fig. 2.8 as a 
function of geocentric distance R for X = 4.65 (i.e., H at 1500 °K) (After F. S. 
Johnson [83]). 




Fig. 2.10. Distribution of H and He in the terrestrial exosphere at 1700 K showing 
the departure from a hydrostatic distribution for H due to escape (after M. Nicolet 

[41]). 
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Fig. 2.11. Fraction of exospheric particles (£) whose velocity v exceeds the escape 
velocity Uoo, as function of the escape parameter X. The critical value of A = 1.5 
represents Opik’s condition for atmospheric “blow-off” [73]. 




Fig. 2.12. Escape parameter X for hydrogen atoms as a function of orbital dis- 
tance for a migrating Jupiter-class (solid line), a Uranus-class (dashed line) and 
an exoplanet with the size of two Earth radii and six Earth masses (dashed-dotted 
line) [84] at XUV energy flux values for G-type stars 1 Gyr after their arrival at 
the ZAMS, where the XUV flux is about 6 times larger than for a 4.5 Gyr old G 
star. One can see that migrating hydrogen-rich exoplanets can easily reach small 
escape parameters due to XUV heating so that exospheric blow-off conditions will 
develop. 
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Fig. 2.13. Escape parameter X for hydrogen atoms for an Uranus-like planet as a 
function of exobase levels in planetary radii (Ro) for different exospheric tempera- 
tures. 



X for hydrogen atoms as a function of orbital distance for migrating exo- 
planets with various masses and sizes. Figure 2.13 shows X for hydrogen 
atoms for an Uranus-like planet as a function of exobase levels in planetary 
radii (Ro) for various exospheric temperatures. Since for exospheric temper- 
atures Too > 10,000 K the exobase level must expand to greater altitudes 
as shown in Fig. 2.14 (cf. Sect. II. 3) the corresponding escape parameter X 
will decrease, reaching the critical value X = 1.5 at altitudes much closer 
to the planet. As the result of escape, the Maxwellian velocity distribution 
is truncated. Thus, at high altitudes the temperature at the exobase is 
not representative. However, a kinetic temperature can be defined by the 
second moment of the velocity distribution function in the moving frame of 
reference, i.e., moving with velocity (vr), according to 

9 T 

T kin = ^J 1(v 2 )-{vr) 2 ]- (2-47) 

The ratio of the kinetic temperature to the temperature at the exobase for 
X c 3> 1.5 is given by [76] 



Tkin _ 7Tl(3 + 2X e )X^ 2 

16V? 5 



(2.47a) 



When the mean energy is dominated by the escaping component, the ratio 
approaches [76] 




(2.47b) 
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100 1000 100C0 tooooo 

Temperature, K 

Fig. 2 . 14 . Rise of exobase levels for an Earth-like planet with present Earth’s 
atmospheric species as a function of exospheric temperature for various XUV energy 
flux values. 1 XUV corresponds to the XUV energy flux of a 4.5 Gyr old solar-like 
G star. One should note that the critical temperature for “blow off” of atomic 
hydrogen on Earth is about 5000 K (courtesy of Yu. N. Kulikov / Polar Geophysical 
Institute, Russian Academy of Sciences, Murmansk, Russia). 



The density distribution of (2.38) applies to a stationary (non-rotating) 
planet. A generalization of the velocity distribution for a corotating exobase, 
satisfying the collisionless Boltzmann equation, leads to a density distribu- 
tion in which the density at the equator exceeds that at the pole, for a 
given planetocentric distance [80]. Further refinements [81], including a non- 
homogeneous (in temperature and density) exobase, lead to even more drastic 
effects of rotation. Features of a rotating versus a stationary exosphere are 
illustrated in Fig. 2.15. The effect of rotation becomes pronounced when the 
rotational parameter of (2.36) Y = QR/uo > 1, i.e., when the rotational 
velocity exceeds the most probable thermal velocity of the exospheric con- 
stituent. 

Table 2.7 lists the values of the rotational parameter at the exobase Y c for 
important constituents in the exospheres of Venus, Earth, Mars and Jupiter. 
Since Y c oc Uq 1 , rotational effects for a given planet are most pronounced 
at low temperatures and for high mass numbers. The upper limits of Y c in 
Table 2.7 are based on estimated exospheric temperatures for minimum solar 
activity. 





II.4 Hydrodynamic Escape and Atmospheric Expansion 



63 




Fig. 2.15. Normalized density distribution for an exospheric constituent ( X = 4.65, 
Y = 1.0) in a rotating exosphere as a function of latitude when the exobase is 
uniform (after R. E. Hartle [81]). 



Table 2.7. Rotational escape parameter. 
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II. 4 Hydrodynamic Escape and Atmospheric Expansion 

If the escape parameter reaches values 1.5 < X < 3 for the main atmospheric 
constituent at the exobase the atmospheric gas starts to expand, so that the 
escape is a process similar in many respects to the supersonic flow of plasma 
from the solar wind [72, 85]. By considering a dynamically expanding non- 
viscous gas of constant molecular weight, in which the pressure is isotropic, 
the steady state equations and energy conservation can be written as 



V • (rijv) — 0, 



(2.48) 



rrijUj (u • VT) + V (njfcpT) = rijmjg 



(2.49) 
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V • (KjVT) = V • 



-rijkBTvj + rij 



mjVj 



— rijmjg • Vj — Qx uv 



(2.50) 



where vj is the bulk velocity of the gas j, g the acceleration of gravity, Kj the 
thermal conductivity parameter and Qxuv is the volume-heating rate caused 
by the XUV radiation. To apply these equations to a dense thermosphere one 
can assume that at a distance 



GMnrij 

UbToXo 



(2.51) 



in the upper atmosphere where no escape takes place (Xq « 30), the temper- 
ature is fixed at T 0 « T e ff and that the atmosphere is sufficiently dense that 
the optical depth of the lower boundary to the XUV flux is much greater than 
1 [85]. Moreover, one can also assume that the XUV flux is absorbed in a dis- 
tance at a somewhat higher altitude rq. Since, there is little energy deposited 
above rq and the gas is continually expanding, the temperature declines as 
r oo. The pressure of the gas also declines toward zero because the gas 
must expand according to a critical solution of the hydrodynamic equations. 
For the estimation of the atmospheric expansion radius rq and the corre- 
sponding atmospheric loss rates one can define the following dimensionless 
parameters [72, 85] for the flux (p 



‘P = L Hj 



k 2 B T 0 

k,q GMrrij 



and for the energy e 



e = ^xuv 



GMrrij 

^bTq kq 



(2.52) 



(2.53) 



where Lb, is the hydrodynamic energy-limited loss rate of constituents j 



_ 47rrj J X uv 
Hj “ GMrm 



(2.54) 



where Ko is a thermal conductivity parameter of ss 4.45 4 erg cm -1 s -1 K -1 
[85]. The corresponding equations for the maximum, escape flux ip n)i and 
associated escape parameter at the atmospheric level r\ where the escape 
takes place, are [72, 85] 






Si + 1 



Xi \( 5 * + 1 )/ 2 






+ 1 



2 



Xo-X! 



(2.55) 



where Xq and X\ are the escape parameters at r — ro and r — r\. One gets 
for Xi 
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Table 2.8. Atmospheric expansion rq for HD209458 b, Jupiter and Saturn as 
function of orbital distance T>, To and corresponding XUV energy flux /xuv with 
Xq = 30. 
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Xi = 




2 

(1 “1“ Tm 




(2.56) 



The simultaneous solutions of these equations for ip m and X\ can be accom- 
plished by numerical methods with specified values of e and Xq. 

Table 2.8 lists the calculated atmospheric expansion ?q in planetary radii 
r p i as function of orbital distance d, temperature To ~ T e ff and corresponding 
/xuv energy flux for a 4.5 Gyr old G-star (the present Sun) for the giant 
exoplanet HD209458 b (Rq = 1.43 ±0.04 Rj up , M 0 = 0.69 ±0.02 Mj up ); [86] 
and for “hypothetical” exoplanets with the exact size and mass of Jupiter 
and Saturn. It should be noted that short-periodic evaporating gas giants at 
very close orbital distances may expand their upper atmospheres up to the 
Roche radius (Roche lobe). In cases where the Roche lobe is lower than the 
hydrodynamic driven expansion ri, the Roche radius limits the expansion 
and the escape leads to a strong mass transfer from the planet to the host 
star [87]. 

One can also see from Table 2.8 that an increase of To tends to decrease 
the expansion. The expansion radii rq decrease if To < T e ff increases, because 
one can see from 2.50 that ro decreases if To grows and Xq is kept as a 
constant value of 30. 

By knowing the expansion radii ?q = (Xo/Xi)ro of the upper atmosphere 
we can calculate the maximum XUV-driven energy- limited loss rate L^. with 
(2.53). After planets have lost their dense hydrogen atmospheres due to the 
energy-limited loss, heavier particles like O, N, O2, N2, CO2, CH2, CH4, NH3 
and other species, which may be present in deeper atmospheric levels or may 
dominate the lower thermosphere, will decrease the energy-limited escape 
rates, since the light hydrogen must diffuse through the heavier constituents 
[88, 89]. 

Analytical escape approximations by planetary winds of atmospheric con- 
stituents can be derived by assuming that the lowermost subsonic regions of 
the hydrodynamic flow controls the differential escape. One should note that 
a real hydrodynamic wind would be transonic and a subsonic approxima- 
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tion is not obviously reasonable, but an analogous treatment for diffusion 
of a heavy trace constituent in a subsonic escaping hydrogen wind provides 
a good approximation to escape rates derived from corresponding transonic 
winds [88]. 

The escape flux 4>j of a heavy atmospheric constituent in a lighter escaping 
atmospheric species by considering ( rrij > rrii) can be written as [88] 



GM ( rrij —rrii) bij 

rgfcsT (! + /,•)’ 



(2.57) 



where fj is the mixing ratio of the heavy gas in the lighter one. It is as- 
sumed that fj is constant over the distances of the density scale height. <fi is 
the escape flux of the lighter atmospheric constituent and bij is the binary 
diffusion coefficient [43] 



k B T 

ki,j 



(2.58) 



with Hij the reduced mass and k{ j the collision coefficient. The escape flux of 
the heavy constituent can also be expressed in terms of the escape parameter 

x, [ 88 ] 



GM (rrij — rrii) bij 
<pir%k B T (1 + fj) 



(2.59) 



or 



g (rrij - rrii) b itj 
4>ikBT (1 + fj) 



The escape parameter Xj cannot be negative and if fj — 0 the heavy com- 
ponent does not escape, then Xj = 0 and 2.59 reduces to 



GM (rrij — rrii) bi 3 __ 

' i= rfMTU + U)' =4>dm 



(2.60) 



so that the escape flux of the lighter constituent in (2.60) is identical to the 
diffusion limited flux 0diff [88, 89]. If the mixing ratio fj <C 1, so that the 
second constituent is no longer major, the presumption that Xj = 0 breaks 
down and the mixing ratio fj must be calculated indirectly for an isothermal 
wind. 

From the equations above one can see that diffusion limited escape is 
generally the case, in which a light atmospheric constituent like hydrogen es- 
capes at the maximum rate for which a heavy major atmospheric constituent 
is retained. Numerical solutions for an isothermal hydrogen wind escaping 
transonically through the hydrostatic background gas must be applied. At 
low temperatures the presence of a second major atmospheric constituent 
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tends to inhibit escape of a lighter trace gas, but under certain conditions its 
escape flux can be enhanced by buoyancy with respect to the heavy major 
species [21, 88, 89]. For a given temperature there exist a threshold mixing 
ratio of the heavy major constituent below which a static solution cannot be 
obtained. Under these conditions diffusion limited flux cannot hold and the 
heavy major atmospheric constituent must also escape. 

Mass fractionation by XUV-driven hydrodynamic escape of hydrogen on 
planets is a promising mechanism for explaining the observed elemental and 
isotopic abundance patterns in terrestrial planetary atmospheres [20, 21, 88]. 
Due to high XUV fluxes of young or active stars (cf. Chap. I) this thermal 
loss process can even remove primitive water oceans from terrestrial planets 

[89]. 

By assuming the amounts of oxygen that are believed to be involved in 
crustal iron oxidation on Venus and Mars, as well as in the massive Venusian 
atmosphere, one finds that primitive water oceans equivalent to, respectively 
0.45-0.2 times the present terrestrial ocean with equivalent depths of 1300- 
600 m could be lost with respectively 30-50% of oxygen initially contained 
in the ocean released to space by XUV-driven hydrodynamic escape [89]. 

Giant gaseous exoplanets with small initial masses close to their host 
stars may even evolve into volatile-rich large terrestrial planets [84] during 
high XUV active periods of their host stars. The remnant cores of giant 
exoplanets orbiting close to their host stars may be detectable for the first 
time with the new generation of space observatories like CoRoT or Kepler. 
Moreover, these considerations show that XUV-driven hydrodynamic escape 
at close orbital distance is a major loss process that should be included in 
any realistic model of the mass and radius evolution of close-in exoplanets. 



II. 5 Atmospheric Sputtering 

Atmospheric sputtering [90, 91] is a collision based energy transfer process 
similar to collisional sputtering from a solid surface [92, 93], which employs 
a so-called collision cascade model. This model is most relevant for cases in 
which the incident ion energy losses by elastic collisions with target nuclei are 
much greater than losses in elastic collisions, which may result in ionization 
or excitation of target particles. In order for solid particles to escape they 
must overcome the surface binding. 

For the escape of atmospheric particles, they must overcome the gravita- 
tional attraction of the atmosphere, where their escape will be unhindered 
by energy degradation through collisions with neighboring particles in the 
exosphere and its lower boundary the exobase. These atmospheric altitude 
regions are equivalent to the first monolayer of material in the sputtering 
of condensed materials. Sputtering from the exosphere is calculated in an 
exact analogy with sputtering from a solid, planetary surface, producing par- 




68 



II Neutral Atmospheres 



Table 2.9. Experimental values of a as a function of mass fraction between the 
sputtered particle and the sputter agent [91, 95]. 
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m B /m A 


0.1 


0.12 
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0.17 
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0.2 
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0.35 


10 


0.5 



tides, which escape directly, or lacking the velocity to escape which will enter 
ballistic trajectories and form an atmospheric corona [91, 94]. 

The sputtering yield Y can be defined as the number of atmospheric 
species ejected per ion incident on an exobase of a gravitationally bound gas. 
This yield includes particles ejected directly in a collision with the ion and 
those ejected due to a cascade of collisions initiated by the ion [94] 



Y = 



0-5<7( Ut > Uoo ) 3dS n 

aj^cosO 7r 2 (cos0) 1 - 6 ?7 oo crD 



(2.61) 



where 6 is the pitch angle of the incident particle, Uoc is the escape energy of 
the particle at the exobase, 0u T >Uoo the c °Ui s i° n cross-section for a particle 
receiving an energy transfer and ctd, is the cross-section for escape of the 
struck particle in order of about 10” 15 cm 2 , according to 

^ < 262 > 



with n c the number density at the exobase, and H the scale height of the 
dominant species at this altitude level. The constants 3/ir 2 and a are obtained 
from the transport equation [94]. One can see from Table 2.9 that the factor 
a depends on the mass fraction mp / m A , where rap is the target particle mass 
and tua is the mass of the incident particle [95] . S n is the so-called stopping 
cross-section [91, 94] 



A 2 

S n = 2? r— - [2 es n (e)} (2.63) 

7^A 

where s(e) is the reduced stopping cross-section with [2ss n (s)] = ln(e) for 
e> 30 and ln( 1 + 1.138e)/(l + 0.0132£-°- 787 + 0.196s 0 5 ) for e < 30. The 
collision cross-section a for a particle receiving an energy transfer Ut > 
is [91, 94] 



da _ 'Ka^f(t 1 / 2 ) 
d t ~ 2t 3 / 2 



(2.64) 
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with the scaling quantities of length au, energy e, A and energy transfer t 



_ 0.8853a 0 

au “ (zpTzp) 



(ifE a) ay 
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Z a ^b^ 2 



g 2 E/ T 

7^A 



(2.65) 

( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 



where ao is the Bohr radius of hydrogen, Z A and Zb the nuclear charge of the 
incident and target particles, 7 the ratio of masses (4m A mB)/(m A + ras) 2 , 
e the elementary charge, £q the dielectricity constant, Ea the energy of the 
incident particle and Ut the energy transfer to the target particle. Since t is 
a function of e, fit 1 ? 2 ) can be evaluated from 



/(<0 = ^ [-n(e)] • (2.69) 

The stopping cross-section S n is between the energy range 0.5-3 keV for 
incident solar wind protons which hit nitrogen or oxygen molecules in the 
order of about 8 x 10 -15 eV cm 2 . The energy transfer cross-section o of 
incident protons with molecular nitrogen or oxygen target molecules is, for 
0.1 keV, about 10 “ 14 cm 2 , for 1 keV about 4 x 10 ~ 15 cm 2 and for 3 keV 
protons about 2 x 10 ~ 15 cm 2 . The escape rate for atmospheric sputtering 5 sp 
from a planetary exosphere who is not protected by an intrinsic magnetic 
field if the incident particles are ions is obtained by [91, 94] 



^sp 



4>i ( 27 tR 2 c ) 



0.5cr(UT> Uce ) 3A,S', 1 

C TdCOsO 7T 2 (cOS0) 1 - 6 J7 oo (7d 



(2.70) 



with R c the exobase distance, and 4>i the incident ion flux which can be solar 
wind protons or back-scattered heavy ions that originate from a planetary 
corona or magnetospheric plasma in the case of a satellite which orbits inside 
the magnetosphere of a giant planet. 



II. 6 Physical Properties of Planetary Atmospheres 

It is commonly accepted that the atmospheres of the terrestrial planets (Mer- 
cury, Venus, Earth and Mars) are of secondary origin, having lost their pri- 
mordial constituents a long time ago. This, however, is not the case for the gas 
giants (Jupiter, Saturn, Uranus and Neptune). The secondary atmospheres 
are most likely due to outgassing of volatile materials from the planetary in- 
terior and radioactive decay products such as He 4 (from the decay of uranium 
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and thorium) and Ar (from potassium K 40 ), with the exception of oxygen in 
the Earth’s atmosphere, which is related to the evolution and presence of life. 
(See Table 2.10 for constituents of planetary atmospheres.) 

The Terrestrial Planets: Mercury, Venus, Earth, and Mars 

Before the Mariner 10 fly-by, it was generally thought that Mercury possessed 
a thin atmosphere similar in composition to the atmospheres of Venus and 
Mars. It was argued that the terrestrial planets should have evolved along 
parallel lines and that a CO 2 atmosphere could be preserved by freezing out 
on the cool darkside. Mariner 10 was therefore designed to search for H, O, C, 
CO and CO 2 , as well as the noble gases He, Ar and Ne. By scaling terrestrial 
noble gas production and diffusion rates and assuming a Lunar-like solar 
wind surface interaction the above mentioned noble gases were thought to be 
the most abundant species [96]. 

The Mariner 10 spacecraft encountered Mercury on 29 March 1974 and 
observed only H, He and O from the UVS instrument [97]. Mariner 10 dis- 
covered also that a direct solar wind interaction was generally not possible, 
since the planet was found to possess a weak magnetosphere [98]. 

A decade after the Mariner 10 fly-bys ground-based observations discov- 
ered the spectrum of Mercury strong emission features at the Fraunhofer 
Na D lines attributed to resonant scattering of sunlight from Na vapor in 
the planetary exosphere [99]. A Na column density of 10 11 -10 12 cm -2 was 
estimated, which corresponds to a Na atom surface density of 2.0 x 10 4 to 
1.5 x 10 5 cm -3 . Since the column content of the observed constituents is of 
the order of about 10 12 cm -2 which is less than the defined column content of 
an exosphere of about 10 14 cm -2 the gaseous envelope of Mercury represents 
an exosphere [100] where the planet’s surface is the exobase. 

Observations of the broad wings in the observed Na D line profile, how- 
ever, suggest that Na extends up to about 700 km above Mercury’s surface 
[99, 101]. This observation indicates the presence of hot Na atoms to match 
the broad wings of the observed line profile. Neither photon-stimulated des- 
orption [102] nor thermal vaporization produce a non-thermal velocity dis- 
tribution enabling Na particles to reach such altitudes. 

Magnetospheric models have shown that Mercury’s magnetosphere could 
be open to the solar wind over substantial areas when the interplanetary 
magnetic field turns southward [103]. Solar wind particles can therefore reach 
the surface where they act as sputtering agents for surface particles. It can be 
shown that particle sputtering could be a source of hot Na atoms with ejection 
speeds larger than 2 km s -1 [94]. Micrometeoroid impact vaporization could 
also be responsible for hot components in Mercury’s exosphere [100, 102, 104, 
105, 106]. 

The Na and O column densities were estimated in earlier work [97, 99, 
104, 107] on the basis of scale heights of about 50 and 80 km, respectively, 
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corresponding to subsolar dayside surface exospheric temperatures of 500- 
575 K. Small concentrations of K in Mercury’s exosphere were also discovered 
[108] so that one can suggest that solar wind ion sputtering might account for 
their existence. Recently Ca was discovered as the sixth element in Mercury’s 
exosphere [109]. 

A more astounding and unexpected discovery was made by the study of 
radar-bright regions near Mercury’s poles, since these areas were attributed 
to H 2 O ice. The radar-bright regions at the poles of Mercury have been used 
to argue for the presence of condensed volatile species in a manner analogous 
to that suggested for permanently shadowed regions on the Moon, although 
it has been suggested that cold silicate minerals would also produce similar 
results. 

In spite of their widely different atmospheres at the present time, it has 
been proposed that the main outgassing constituents from the planetary in- 
teriors have been the same for Venus, Earth and Mars, i.e., water vapor and 
CO 2 with a small amount (<1%) of N 2 . Their different evolutionary paths 
have been suggested to be the consequence of the initial effective surface 
temperature due to their distance from the Sun together with a “runaway” 
greenhouse effect on Venus due to H 2 O and CO 2 [110, 111, 112]. Venus has 
a dense and hot atmosphere, while Mars has a thin and cold one relative 
to Earth [113]. The surface pressures of Venus, Earth and Mars are in the 
approximate ratio 100 : 1 : 0.01, while their dayside surface temperatures T s 
are about 750 K, 300 K and 250 K, respectively. Both Venus and Mars have 
CO 2 as their principal atmospheric constituent, whereas Earth has N 2 and 
O 2 [48]. The temperature and pressure distribution for Venus, Earth and 
Mars are shown in Figs. 2.16 and 2.17. It has been suggested that because 
of the initial high surface temperature of Venus, H 2 O, if present at the for- 
mation of the planet, was always in the atmosphere in the form of vapor, 
since the surface temperature was always above the boiling point of water at 
the given pressures in the early history of the planet. The complete absence 
of liquid water together with the high surface temperature then allowed a 
substantial accumulation of CO 2 in the atmosphere, since reactions between 
silicates in the crust and CO 2 proceed rapidly only in the presence of liquid 
water and are also strongly temperature dependent. The present lack of H 2 O 
is explained as either due to photodissociation by solar ultraviolet radiation, 
with the hydrogen having escaped and the oxygen consumed by various ox- 
idation processes at the surface, having escaped or lacking due to an initial 
absence of H 2 O. 

The composition of the Earth’s atmosphere is thought to be due to the 
initial (lower) surface temperature at an orbit inside the liquid water habit- 
able zone, allowing water to condense at the surface and accelerating chemical 
weathering and carbonate formation of atmospheric CO 2 and the accumu- 
lation of the inert gas, nitrogen, supplied from volcanic activity and oxygen 
associated with life. 
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Fig. 2.16. Temperature distributions for Venus, Earth and Mars consistent with 
our present understanding of these planetary atmospheres. For Venus the low night- 
side exosphere temperature results from the slow retrograde rotation rate. 




Pressure (atmospheres) 



Fig. 2.17. Model pressure distributions for the atmospheres of Venus, Earth and 
Mars based on observations and theory. 



Mars, on the other hand, with its low initial temperature, caused the vol- 
canic steam to freeze at the surface and allowed the accumulation of CO 2 
in the atmosphere, but to a much lesser extent on account of its low tem- 
perature [113]. Although the above evolutionary path of the atmospheres of 
the terrestrial planets is largely speculative, it is reasonably consistent with 
present ground based and spacecraft observations. 
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The physical parameters of Earth have, of course, been well established by 
observations. In addition to the principal atmospheric constituents and their 
dissociation products, the upper atmosphere of Earth contains helium (which 
is of radiogenic origin) and hydrogen (which is the result of photodissocia- 
tion of water vapor and methane) as the major light constituents, in addition 
to atomic oxygen. Because of the range of exospheric temperatures over the 
solar cycle, the terrestrial hydrogen corona is variable. A hydrogen corona 
has also been observed for Mars [113] and Venus [111, 112]. Helium, though 
probably present, particularly on Venus, has yet to be identified directly. The 
presence of oxygen in the upper atmospheres of both the Venus and Mars 
has been established, but its concentration is extremely small compared to 
Earth, raising the interesting question of the stability of a CO2 atmosphere. 
Since photodissociation should be very effective, competing processes of re- 
combination exist if the main heavy constituent of the upper atmospheres of 
Mars and Venus is indeed CO2. The CO2 problem remains one of the most 
enigmatic ones regarding the atmospheres of Mars and Venus [114]. 

Because of the relatively small scale heights the exobase for Mars and 
Venus occurs at a much lower altitude (~ 220 km) than for Earth (~ 500 km). 
Due to the absence of a significant planetary magnetic held, the solar wind 
can therefore interact with the exospheres of Venus and Mars at rather low 
altitudes. 



The Gas Giants: Jupiter, Saturn, Uranus, and Neptune 

In contrast to the terrestrial planets, the gaseous giants or outer planets 
have reducing atmospheres with large amounts of hydrogen, suggesting that 
they represent primordial atmospheres whose origin is similar to that of the 
Solar System. According to spacecraft observations, Jupiter [115] and Saturn 
[116] have atmospheres which in addition to hydrogen contain helium (He), 
methane (CH4) and ammonia (NH3), while Uranus [117] and Neptune [118] 
are richer in the heavier elements, which may have their origin from their 
icy-rocky cores. Although, the composition of the material in the deep layers 
of Uranus and Neptune is not well known, it must be much denser than 
hydrogen or He, but not as dense as Mg-Si and rocky iron material which 
compromises terrestrial planets that formed in the inner Solar System. It is 
commonly believed that this material may largely consist of abundant ices 
which contain H2, CH4 and NH3 and rock, which is a solar mix of those 
elements, which are the major constituents of rock on Earth [119, 120, 121]. 
The compositions of this rock may be 38% SiCU, 25% FeS and 12% FeO [120] 
or 39% Si0 2 , 32% Fe, 27% MgO and 2% Ni [119, 120]. 

As the result of their great distance from the Sun, the exospheric tempera- 
ture contribution from the XUV radiation of the outer planets is close to their 
effective temperature (~ 125 K for Jupiter) [70], but there are other sources 
of thermospheric heating such as dissipation of acoustic energy, particle heat- 
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ing, and atmospheric gravity waves, which result in exospheric temperatures 
between 700 and 1000 K. 

The Reduced Atmospheres of Titan, Triton, and Pluto 

Saturn’s and Neptune’s large satellites Titan and Triton as well as the planet 
Pluto have reduced atmospheres consisting of CH 4 and N 2 and at confirmed 
in Titan’s case organic molecular compounds like HCN. Generally these small 
bodies may represent small cores of Uranus-type planets. Bodies like Titan 
are also thought to have a rocky core, which is hidden below a H 2 O ice layer 
with trapped CH 4 , NH 3 and other volatiles. Recent ground based observa- 
tions with the United Kingdom Infrared Telescope (UKIRT) support this 
suggestion, because the characteristic of H 2 O ice was discovered on Titan’s 
surface in narrow wavelengths centered at 0.83, 0.94, 1.07, 1.28, 1.58, 2.0, 2.9 
and 5.0 micrometers [122]. During the formation of Saturn and Titan the H 2 O 
ice-clathrate melted, so that the trapped CH 4 and NH 3 could be released in 
large amounts and built up the atmosphere. The outgassed CH 4 reacted with 
NH 3 and produced Titan’s dense N 2 /CH 4 atmosphere. After the body cooled 
out, the orbital distance of about 10 AU is responsible for the present surface 
temperature of about 100 K so that the H 2 O is in a frozen state. On the 
other hand CH 4 and C 2 H 6 may be liquid on the surface and thin haze layers 
may cover Titan’s icy surface and influence the albedo at short wavelengths, 
while leaving its 2.9 and 5.0 micrometer reflectivities unaffected [122]. 

Figure 2.18 shows the pressure and temperature distribution in the atmo- 
spheres of Venus, Earth, Mars, Jupiter, Saturn, Uranus, Neptune, Pluto and 
Saturn’s large satellite Titan. 



Table 2.10. Main gases identified in planetary atmospheres. 



Mercury Na, O, K, Ca, H, He, ? 

Venus C0 2 , N 2 , S0 2 , H 2 S0 4 , CO, H 2 0, O, H 2 , H, D 

Earth N 2 2 , 0 2 , H 2 0, Ar, C0 2 , Ne, He, CH 4 , K, N 2 0, H 2 , H, O, 0 3 , Xe 

Mars C0 2 , N 2 , 0 2 , CO, H 2 0, O, He, H 2 , H, D, O 3 

Jupiter H 2 , He, H, CH 4 , NH 3 , CH 3 D, PH 3 , HD, H 2 0 

Saturn H 2 , He, CH 4 , NH 3 , CH 3 D, C 2 H 2 , C 2 H 6 

Uranus H 2 , He, CH 4 , NH 3 , CH 3 D, C 2 H 2 , 

Neptune H 2 , He, CH 4 , NH 3 , CH 3 D, C 2 H 2 , C 2 H 6 , CO 
Pluto N 2 , CH 4 , ? 

Titan N 2 ,CH 4 , HCN, organics 
Triton N 2 , CH 4 , ? 
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Fig. 2.18. Temperature versus pressure distribution in the atmospheres of the 
planets in the Solar System. 



Exoplanets 

The observation of the first giant exoplanet surrounding a solar-type main 
sequence star was announced in 1995, based on very precise radial velocity 
(Doppler) measurements. Before 1995 the question whether planets existed 
outside our solar system was a hypothetical one. Today (February 2004), we 
know 119 planets in 104 extrasolar systems, where 13 are multiple systems 
and which are investigated by various scientific groups all over the world 
[123]. With the discovery of planets outside our own Solar System modern 
science is continuing the Copernican revolution, placing our planetary system 
among a wealth of systems existing in the Universe. The exoplanetary systems 
discovered so far show a surprising diversity of orbital parameters. Most of 
these systems do not resemble our own Solar System. There is presently no 
consensus as to why these planetary systems are so diverse, but it seems that 
migration due to interactions between planets and the protoplanetary disk 
are an important ingredient [124, 125]. 

Figure 2.19 shows a comparison of 119 giant exoplanets as a function of 
Jovian mass and distance with the planets of our Solar System. Bodies with 
masses around and above 15 Jovian masses are brown dwarfs where fusion 
of deuterium isotopes can occur. The difference between a brown dwarf and 
a low mass M star is that the mass of brown dwarfs is so low that hydrogen 
fusion does not start. Among the large number of giant exoplanets detected 
up to now there is a subset of exoplanets in very close orbits around their host 
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D [AU] 

Fig. 2.19. Comparison of the discovered exoplanets in Jovian mass and orbital 
distance with Mercury, Venus, Earth, Jupiter, Saturn, Uranus and Neptune. The 
dotted line shows the orbital distance where giant hydrogen-rich exoplanets can 
experience hydrodynamic conditions. The dashed line marks the border where the 
mass loss of gas giants can be very effective. 



stars. Through transit observations it was shown that at least two of them, 
namely HD209458 b and OGLE-TR-56 b, orbiting their hosts stars at 0.045 
AU and 0.0225 AU, respectively, are known to be Jupiter-like gas giants with 
hydrogen as their main atmospheric constituent. The observed parameters 
for HD209458 b are R 0 = 1.43 ± 0.04 R Jup , M 0 = 0.69 =b 0.02 M Jup , r = 5.2 
Gyr [86] and for OGLE-TR-56 b, R 0 = 1.3 ± 0.15 R Jup , M 0 « 1.6 M Jup , 
and r = 3 ± 1 Gyr [126, 127]. For all the other detected exoplanets only an 
estimated mass is known due to radial velocity measurements. 

The recent observation of an extended and hydrogen-evaporating atmo- 
sphere for HD209458 b due to absorbtion in the Lyman-a line with the HST 
raises the question of atmospheric stability against escape at very short or- 
bital distances [72, 128]. 

By solving (2.55) and (2.56) for getting the atmospheric expansion radius 
ri one obtains for HD209458 b a maximal energy-limited hydrodynamic loss 
rate of « 10 12 gs _1 , which is in good agreement with the value determined 
from planetary transit observations [128]. Indeed, a lower limit to the escape 
rate of « 10 10 g s _1 that can be orders of magnitude higher because of the 
saturation of the absorption line [72] was estimated. The dotted line Fig. 2.19 
shows the border where giant hydrogen-rich exoplanets begin to experience 
XUV-driven hydrodynamic expansion and high loss rates. At orbits closer 
than 0.07 AU (dashed line) the hydrogen loss rates can be very high so that 
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Table 2.11. Main gases identified in the atmospheres of brown dwarfs and of the 
short-periodic giant exoplanet HD209458 b. 

Brown dwarfs CO, CH 4 , NH 3 , Na, H 2 , H, ? 

HD209458 b (“Hot Jupiter”) H 2 , H, Na, C, O, ? 



exoplanets may evaporate a large fraction of their atmospheric mass, even 
down to their core size during the early XUV active period of their host 
stars. 

By comparing the spectra taken during the transits of the atmosphere of 
HD209458 b by several transits with the Space Telescope Imaging Spectro- 
graph (STIS) instrument on the HST, with those taken at other times, it was 
found that traces of atomic Na in the spectra of the transits could have their 
origin only in the atmosphere of the exoplanet, although Na is not the major 
atmospheric component [129]. After the observation of a hydrodynamically 
driven strong planetary hydrogen wind [128] it is clear that the main atmo- 
spheric constituent of HD209458 b will be like that on the gas giants of our 
Solar System, H 2 . Recent observations with the STIS of transits of HD209458 
b revealed also the loss of heavy species C and O from the upper atmosphere 
of HD209458 b [130] as shown in Table 2.11. These heavy species, however, 
must be carried away from the planet due to the strong hydrodynamic-driven 
hydrogen wind as discussed in this chapter. For comparison, the escape rate 
resulting from the Jeans formula at a temperature equal to T e // is < 1 g s -1 . 
For the closer exoplanet OGLE-TR-56 b, with a measured radius, one can es- 
timate a present mass loss rate of « 5 x 10 12 g s -1 . Due to much higher XUV 
fluxes after their host stars arrived at the ZAMS both exoplanets may have 
lost by evaporation a significant fraction, up to twice their original masses, 
during their lifetime [72]. 

The first transit observed (HD209458 b) was not a detection, since the 
planet was previously found by radial velocity searches. The OGLE project 
may have detected the first planet by the transit method (OGLE-TR-56 b). 
The reason for the small detection rates so far are mainly observational con- 
straints. These will be overcome with the near- future space missions C 0 R 0 T 
(ONES), Kepler (NASA) and Eddington (ESA), which are expected to detect 
many gas giants and small terrestrial-like exoplanets by the transit method. 

It also seems reasonable to assume that planets resembling Uranus and 
Neptune, or slightly less massive ones, may have formed in cold regions of 
a protoplanetary disk and migrated inward, possibly into the habitable zone 
of the host star where liquid water can be present at their surface. These 
planets, termed ocean planets [84], would be more interesting as their large 
radius makes them more easily detectable by transit missions like C 0 R 0 T, 
Kepler and Eddington. 

On such ocean planets, expected primary volatiles would be H 2 O, NH 3 
and CO 2 , because the assumed initial composition of ice is similar to that of 




78 



II Neutral Atmospheres 



comets, e.g. w 90% H2O, 5% NH3 and 5% CO2 [84]. A two Earth size and 
six Earth mass ocean planet can release hydrogen via photodissociation of 
NH3 to an amount of about 1.6 xlO 23 kg, corresponding to an atmospheric 
pressure of about 1.2 GPa which can be lost in a very short time (^25 Myr) 
by energy-limited hydrodynamic escape during the first Gyr when the planet 
orbits distances < 1 AU from a young solar-like star due to a 100 times 
higher XUV flux during the first 0.1 Gyr after the star arrived at the ZAMS. 
After the majority of the hydrogen is lost, the heavier atmospheric species 
like CO2, O and N become relevant in the upper atmosphere. 

It seems logical that the atmospheres of exoplanets should consist of simi- 
lar constituents to the planets in the Solar System. Future space missions like 
Darwin (ESA) or TPF (NASA) which will have the capability to study the 
atmospheres of exoplanets by thermal emission in the infrared or by reflected 
light in the visible should be prepared to observe hydrogen- dominated pri- 
mordial atmospheres of exoplanet ary gas giants, secondary atmospheres of 
terrestrial planets and reduced atmospheres of migrated more exotic sub- 
Uranus type bodies or atmospheres, which are outgassed by the cores of 
evaporated hydrogen-rich Uranus-class exoplanets. 




Ill Thermal Structure of Upper Atmospheres 



III.l Departures from Thermal Equilibrium 

Electrons released in the photoionization process due to the absorption of 
solar XUV radiation in planetary atmospheres may have initial kinetic ener- 
gies of one 100 eV. The energy E of photoelectrons resulting from photons 
of energy \)v interacting with an atmospheric constituent whose ionization 
potential is IP is given by 



E = f)i/ - IP . 



These photoelectrons lose their excess energy by both inelastic and elastic 
collisions [131]. As long as the energy of the primary photoelectron exceeds 
IP, secondary ionization may occur, while below IP energy loss will occur 
by excitation of atmospheric constituents. The rate of energy loss per unit 
pathlength due to an inelastic collision with a neutral constituent whose 
density is 5 can be expressed by 



J_ d E 
nj dx 



E'oj 



(3.1) 



where E' is the excitation energy and <jj is the appropriate cross-section. 

For the principal atmospheric constituents of the terrestrial planets, pho- 
toelectron energies E > 20 eV lead to optical excitation. For photoelectron 
energies 5 < E < 20 eV, electronic excitation becomes the predominant en- 
ergy loss process. Below 5eV vibrational and, at the lowest energies, rota- 
tional excitation of molecular species leads to photoelectron energy loss by 
inelastic collisions with neutral particles. However, at energies below 2eV, 
elastic Coulomb collisions are most important for the thermalization of pho- 
toelectrons. The rate of energy loss per unit path length for elastic collisions 
is given by 



1 d E 2 m e m 

N dx (m e + m) 2(JD 



(3.2) 



where N is the number density and m the mass of the collision partner 
and ctd is the collision (momentum transfer or diffusion) cross-section [9]. 
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Due to the mass dependence, photoelectrons will lose their energy in elastic 
collisions primarily to the ambient electrons, rather than to the ions, leading 
to a preferential heating of the electrons gas, and an electron temperature 
greater than that of the heavy particles (ions and neutrals). 

This can also be seen by considering the time constant for Maxwellian- 
ization of electrons or ions by self-interaction, according to 



.£§ 

~N~ 



(3.3) 



Numerically, the self-collision time is given by [132] 



11.74m) T§ 

Tjj = NjZ 4 lnA ’ 

where rhj is the mass in amu, and Z the charge of the particles whose density 
is N and temperature T. 

Since ra e = l/1836ra p , the self-collision time for electrons is by a factor 
1/43 less than that for protons; In A is a parameter called the Coulomb log- 
arithm depending on T and N , having numerical values between 10 and 20 
for typical conditions in planetary ionospheres (cf. Chap. V). 

Equipartition between the hotter electrons and the cooler ions will then 
take place by electron-ion Coulomb collisions. 

The equipartition time, r e i is defined according to 



d Te 
d t 



Te-Tj 

^ei 



and is related to the electron self-collision time r ee by 



( 3 . 4 ) 



^ 1 rh[ _ 5.7m e rhi / T e T\ \ 2 
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Since r e i r ee , the electron temperature T e can exceed the ion temperature 

T[ in altitude regions where the main cooling process for electrons is that of 
electron-ion collisions, as long as the electron gas is selectively heated. 

The detailed thermal balance of a planetary ionosphere and the individual 
electron and ion temperature depend on the sources and sinks, both local and 
non-local, of the thermal energy of the ionospheric particles [133, 134, 135]. 
These topics will be discussed in the following sections. 



III. 2 Electron and Ion Temperatures 

The temperature of the electrons and ions of a planetary ionosphere is gov- 
erned by heat balance equations of the form 
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as long as the energy of any bulk motion is small compared to that of the 
random thermal energy. The subscript “e” and “i” refer to the electron and 
ion properties, respectively. K refers to the thermal conductivities and Q and 
L refer to the rates of volume heat production and loss. If the ionosphere is 
pervaded by a magnetic field, the heat conduction term is strongly controlled 
by the direction of this field; the sin 2 / factor, where I is the inclination 
(dip angle) of the magnetic field, implies that heat conduction is important 
primarily along magnetic flux tubes, but is strongly inhibited across the field. 



Electron Heating 

As stated in the previous section, the initial source of energy for heating the 
charged particles in the ionosphere derives from the photoionization process 
in which energetic photoelectrons are created which eventually thermalize. 
This leads to an enhanced electron temperature; the difference between elec- 
tron and ion (and neutral) temperatures, i.e., the absence of thermal equi- 
librium in a planetary ionosphere, is controlled by the detailed balance of 
sources and sinks of thermal energy for both the electrons and ions. The 
initial photoelectron energy spectrum, which depends on the characteristics 
of the ionizing solar XUV radiation and the composition and density distri- 
bution of the neutral atmosphere and appropriate cross-sections, represents 
the starting point for the ionospheric heat input. The local heating rate Q e 
is related to the local ionization rate q by 

Qe = sq ( 3 - 6 ) 



where e is the energy per ion pair formation which is transferred locally to 
the electron gas in form of heat and 



eq = / q{E) dE/[ 1 + (dE/dx)J{dE/dx) e ] 

JE th 

where E is the initial photoelectron energy, E t h is the thermal energy of 
the electron gas and the subscripts “n” and “e” refer to energy loss rate to 
neutrals and electrons, respectively. 

From the previous definition of the ion production rate (cf. Chap. I) it 
can be seen that the local heating rate has an altitude dependence which is 
largely controlled by that of the ion production rate, and is thus given by 



Q(z) = £(z)n 0 <Ti$xuv exp 



z ' 

H. ' 



At altitudes where inelastic collisions of electrons with neutrals are impor- 
tant, the heating efficiency is roughly inversely proportional to the neutral gas 
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density, while at higher altitudes e is strongly dependent on the electron, ion 
and neutral density. At these higher altitudes, specifically, where the electron 
mean free path A e ,n becomes comparable to the scale height of the neutral gas 
responsible for scattering, the photoelectrons can escape, i.e., travel signifi- 
cant distances before becoming thermalized [133]. This photoelectron escape 
flux is therefore responsible for non-local heating of the electron gas. The 
escape level for photoelectrons, in analogy to the definition of the neutral 
exobase, is defined by the condition A e , n = H or n = ( a sc H)~ l where a sc is 
the scattering cross-section for electrons in the neutral gas (<j sc = 10 -15 cm 2 ). 
The magnitude of the photoelectron escape flux can be estimated from simple 
scattering considerations [135]. The continuity equation for the photoelectron 
flux (/) can be expressed by 



^ = (3.7) 

where q is again the production rate of photoelectrons and L represents the 
loss due to electron-neutral particle scattering along the path s. 

The factor 1/2 is based on the assumption of an initial isotropic photo- 
electron distribution, so that one half of the electrons produced at a given 
level where escape can take place are directed upward. If there is a planetary 
magnetic field pervading the ionosphere, the path of the photoelectrons is 
constrained by field lines and the scattering also depends on the local pitch 
angle a. For a probability of 0.5 for an inelastic collision to remove a photo- 
electron from the flux, the loss rate can be expressed by 

L — \n<j sc (j)seca . 

Integration of (3.7) gives for the unscattered photoelectron escape flux for 
the limit (z oo) 



or (for a = 0), 



0oo — Q^e,n COS OL 



4>oc = — £oc = — . (3.8) 

°sc Vsc 

With ionization rate coefficients J = 10 -7 sec -1 and scattering cross- 
sections <r sc = 10“ 15 cm 2 , photoelectron escape fluxes of the order of </> = 
10 8 cm -2 sec -1 are obtained on Earth, and photoelectron energy fluxes of 
the order of 10 9 eV cm -2 sec -1 , assuming an average photoelectron energy of 
10 eV. Such energy fluxes can lead to heating rates at high altitudes which are 
orders of magnitudes larger than those resulting from local ion production. 
Thus, escaping photoelectrons can represent an important source of heating 
in the upper portions of planetary ionospheres, and can also lead to heating 
of nighttime ionosphere, if the magnetically conjugate point is already sunlit 
[136, 137]. The heating rate due to a flux of photoelectrons is given by 
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where d E / dz is the energy loss of photoelectrons due to elastic collisions with 
ambient electrons, leading to heating of the electron gas. The energy loss of 
photoelectrons due to elastic collisions has usually been determined from the 
formula of Butler and Buckingham [131, 138, 139, 140, 141] 
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where 

ff ( u ) = 4 = ( r e “* 2 dx - 2ue ““ 2 

v 71 * \J 0 

v t h is the electron thermal speed and In A is the Coulomb logarithm. For 
energies E > 3 eV this equation can be approximated by [135] 

^ = -^L^e [eVsec- 1 ]. (3.11) 

dt E 2 

In the presence of a magnetic field the right-hand side of (3.11) is modified 
by a factor (cos a sin/) -1 . The Butler and Buckingham formula neglects the 
energy loss of photoelectrons due to the generation of plasma waves by the 
Cerenkov mechanism. The complete expression for the energy loss of fast 
photoelectrons to thermal electrons is given by [138, 139] 
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where &(v/v t h) represents the energy loss due to plasma wave generation and 
<jn is the plasma frequency c^n = (dirNe 2 /m e ) 2 . For energies E fepT, this 
formula can be approximated by 



dE 

dt 



,2 ^2 



In 



• < 



mv 

7e 2 o;N 



for ksT <C E < 



me 



4 ^ 



2h 2 



In 



V^N 



for E > 



me 



(3.13) 



where 2nh is Planck’s constant and 7 = 1.781 is Euler’s constant. Energy 
loss rates based on the Butler and Buckingham formula, as applied by many 
authors, may be as much as 70% lower than those based on (III. 2) [138, 140]. 
A useful analytic form of (III. 2) is given by [141] 
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